
LABORATOIRE DE PHYSIQUE THEORIQUE 
UNITE MIX I E DE RECHERCHE n° 7085 
UNIVERSITE LOUIS PASTEUR 
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 
Institut de Physique - 3 rue de 1'Universite- 67084 Strasbourg Cedex 
Tel : 03 90 24 07 38 / Fax : 03 90 24 06 54 




Effects of a strict site-occupation 
constraint in the description of 
quantum spin systems at finite 

temperature 



These 



pour l'obtention du grade de DOCTEUR DE L'UNIVERSITE LOUIS PASTEUR 

STRASBOURG I 
Discipline : Physique Theorique 



Presentee par 

Raoul Dillenschneider 

Le 8 Septembre 2006 a Strasbourg 



Directeur de these 

Jean Richert, Directeur de Recherche 

Jury 

Rapporteur interne : Daniel Cabra, Professeur 
Rapporteur externe : Frederic Mila, Professeur 
Rapporteur externe : Pierre Pujol, Maitre de conferences 
Examinateur : Malte Henkel, Professeur 

Examinateur : Michel Rausch de Traubenberg, Maitre de conferences 



Abstract 



We study quantum spin systems described by Heisenberg-like models at finite tempera- 
ture with a strict site-occupation constraint imposed by a procedure originally proposed 
by V. N. Popov and S. A. Fedotov [HE!- We show that the strict site-occupation con- 
straint modifies quantitatively the behaviour of physical quantities when compared to the 
case for which this constraint is fixed in the average by means of a Lagrange multiplier 
method. The relevance of the Neel state with the strict site-occupation contraint of the 
spin lattice is studied. With an exact site-occupation the transition temperature of the 
antiferromagnetic Neel and spin liquid order parameters are twice as large as the critical 
temperature one gets with an average Lagrange multiplier method. We consider also a 
mapping of the low-energy spin Hamiltonian into a QED 3 Lagrangian of spinons. In this 
framework we compare the dynamically generated mass to the one obtained by means of 
an average site-occupation constraint. 



Resume 



Nous etudions des systemes de spin quantiques a temperature finie avec une contrainte 
d'occupation stricte des sites au moyen d'une procedure introduite par V. N. Popov 
et S. A. Fedotov [HE]- Nous montrons que cette contrainte modifie le comportement 
d'observables physiques par rapport au cas ou cette contrainte est fixee de fagon moyenne 
par la methode des multiplicateurs de Lagrange. La pertinence de l'etat de Neel est 
etudiee en presence de la contrainte stricte d'occupation des sites du reseau de spin. La 
temperature de transition des parametres d'ordre antiferromagnetique de Neel et d'etat 
de liquide de spins sont doubles par rapport a ceux obtenu par la methode moyenne des 
multiplicateurs de Lagrange. Nous considerons l'Hamiltonien de basse energie decrit par 
un Lagrangien de QED 3 pour les spinons. Dans ce contexte la masse generee dynamique- 
ment est comparee a celle obtenue par la methode d'occupation moyenne de site. 
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Chapter 1 
Resume 



1.1 La procedure de Popov et Fedotov (PFP) 

Apres la decouverte de la supraconductivite par Bednorz et Miiller [TU] un effort theorique 
important fut consacre a la recherche d'une explication a ce phenomene qui se manifeste 
a une temperature critique elevee. De nombreux supraconducteurs furent decouverts 
parmi lesquels on retrouve les cuprates avec La 2 - x Sr x CuO^ et Bi 2 Sr 2 CaCu 2 0s qui sont 
des exemples de supraconducteurs dopes en trous, les ruthenates Sr 2 RuO i) les metaux 
a fermions-lourds tels que UPt 3 et ZrZn 2 , ainsi que les materiaux organiques comme 
k - (BEDT - TTF) 2 - Cu[N(CN) 2 ]Br. 

Tous les cuprates partagent la meme structure atomique composee de couches de 
Cu0 2 , tenues pour responsables de la formation de paires de Cooper et de la supracon- 
ductivite, et intercalees de couches de substance dopante et/ou non-dopante. 

Dans ce manuscrit nous nous concentrerons sur les couches Cu0 2 et plus precisement 
sur la phase isolante antiferromagnetique des supraconducteurs haute temperature. En 
effet, la phase isolante des materiaux comme les cuprates peut etre modelisee par le modele 
de Heisenberg. 

L 'etude de la phase isolante est motivee par le fait que les correlations sous-jacentes a la 
phase antiferromagnetique isolante pourrait se prolonger dans la phase supraconductrice 
sous 1' effet du dopage. En d'autre termes nous attendons que la phase supraconductrice 
garde quelques traits de caractere de la phase antiferromagnetique isolante. 

Notre travail est centre sur l'etude des effet s de la contrainte stricte d'occupation de site 
de spin caracterisant la phase isolante dans la description de systeme de spin quantique a 
temperature finie. Cette contrainte consiste a imposer exactement un spin S = 1/2 pour 
chaque site de reseaux. 

L'implementation d'une telle contrainte fut introduite par Popov et Fedotov [HE] . Con- 
trairement a d'autre methodes basees sur l'utilisation de multiplicateurs de Lagrange 
|H1 I7T| . conduisant a une occupation moyenne, la procedure de Popov et Fedotov 
(PFP) evite ce traitement approximatif au moyen de l'introduction d'un potentiel chim- 
ique imaginaire. Nous portons notre interet a l'analyse et l'application de la PFP sur des 
systemes antiferromagnetiques de spins a temperature finie, et nous la comparons aux 
resultats obtenus par la methode des multiplicateurs de Lagrange. 
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Chapter 1. Resume 



1.2 Proprietes magnetiques du modele de Heisen- 
berg par la PFP 

Le chapitre |U presente et discute l'application du champ moyen et du developpement en 
nombre de boucles pour la determination de proprieties physique de systemes antiferro- 
magnetique de type Heisenberg pour des dimensions D du reseau [2*T] . 

Des travaux recents sur des systemes de spins quantiques discutent l'existence possible 
d'etats de liquide de spin, et pour deux dimensions d'espace de la competition ou transi- 
tion de phase entre etats liquide de spin et antiferromgnetique de Neel [2E1 EH EH EH] . II 
est egalement connu que les supraconducteurs (cuprates) presentent une phase antiferro- 
magnetique [HO] , 

Dans le chapitre 0] nous focalisons notre attention sur la phase dont le champ moyen est 
de type Neel dans la description des systemes quantiques de spin represented par le modele 
de Heisenberg. Plus precisement, nous presentons une etude detaillee de l'aimantation 
et de la susceptibilite magnetique pour ces systemes de reseaux de spin en dimension D 
et a temperature finie. Le but de ce travail est d'etudier la pertinence de V ansatz de 
Neel comme approximation de champ moyen en utilisant la PFP et pour les intervalles 
de temperatures < T < T c ou T c est la temperature critique. Dans l'objectif d'obtenir 
une reponse sur ce point nous calculons les contributions des fluctuations quantiques et 
thermiques au-dela de 1' approximation de champ moyen et sous la contrainte d'occupation 
stricte d'un seul spin par site de reseau [5J I2H 1221 H21 EH] • Une comparaison est effectuee 
entre nos resultats et ceux obtenus par la theorie des ondes de spins sur les modeles de 
Heisenberg et XX Z. 

II est montre que la PFP introduit un grand decalage de la temperature critique par 
rapport a celle obtenue par la methode ordinaire des multiplicateurs de Lagrange. En 
effet, il apparait un doublement de la temperature de transition entre un etat antiferro- 
magnetique de Neel et la phase paramagnetique [3 122] ■ Nous montrons dans le chapitre 
IH et c'est la une des contribution originales de ce travail de these, qu'a basse temperature 
l'aimantation ainsi que la susceptibilite magnetique sont egalles en valeur a celles obtenues 
par la theorie des ondes de spin, comme le montre la figure 14.31 du chapitre 0] La figure 
!4.Hl represente l'aimantation en champ moyen (en traits tiretes), l'aimantation avec fluc- 
tuation en utilisant la PFP (en traits plein) et enfin l'aimantation calculee par la theorie 
des ondes de spins (en traits pointilles), pour des systemes tri-dimensionnels. La superpo- 
sition des courbes a basse temperature apparait clairement sur cette figure. Le decalage 
entre l'aimantation moyenne et celle obtenue en considerant les fluctuations est due a 
l'effet intrinseque des fluctuations quantiques et thermiques sur la moyenne statistique 
des orientations de spin. La PFP n'est pas responsable de ce decalage. 

A plus haute temperature les contribution des fluctuations de nature quantiques et 
thermiques croissent en une singularity au voisinage de la temperature critique. L'hypothese 
que le champ moyen de Neel contribue pour une majeure partie a l'aimantation et a la 
susceptibilite n'est plus valide. En approchant de la temperature critique T c l'aimantation 
en champ moyen tend vers zero et les fluctuations croissent de plus en plus a l'ordre d'une 
boucle. Ceci resulte en une divergence de la valeur theorique prevue pour l'aimantation 
globale du systeme de spins. Ce comportement est commun aux modele de Heisenberg 
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et XX Z. Nous en deduisons que la brisure de symetrie induite par le choix de l'etat de 
Neel n'est pas impliquee dans l'apparition des divergences calculees. 

L'influence des fluctuations decroit avec l'augmentation de la dimension D du systeme 
due a la reduction des fluctuations par rapport au champ moyen. 

En dimension D = 2 l'aimantation verifie effectivement le theoreme de Mermin et 
Wagner Pour T ^ 0, les fluctuations sont plus importantes que la contribution du 
champ moyen. Ainsi dans une description plus realiste de la physique il est necessaire 
de prendre en consideration d'autre champ moyen. En effet, Ghaemi et Senthil [2E] ont 
montre, a l'aide d'un modele specifique, qu'une transition du second ordre d'un etat de 
Neel vers un liquide de spins pourrait apparaitre en fonction de la valeur des couplages. 
Ceci nous conduit au chapitre dans lequel une etude de l'influence de la PFP est menee 
sur differents ansatz de champ moyen et pour des reseaux bi-dimensionnels. 

Le point originale de ce chapitre |U est l'utilisation de la procedure de Popov et Fedotov 
dans la calcul de proprieties magnetique a temperature finie et pour un ordre a une boucle 
en perturbation [21]. 

1.3 Champ moyens pour le modele bidimensionnel 
de Heisenberg 

Nous resumerons dans cette section le contenu du chapitre ou nous considerons la PFP 
appliquee a des systemes de spin quantique a temperature finie. Le but de ce chapitre est 
de confronter l'approche utilisant la PFP a celle utilisant un multiplicateur de Lagrange 
pour fixer la contrainte d'occupation par site du reseau de spin. 

La description des systemes quantiques de spins a temperature finie passe generalement 
par l'utilisation d'une procedure de point-selle qui est une approximation d'ordre zero de 
la function de partition. La solution de champ moyen ainsi generee peut fournir une 
approximation realiste de la solution exacte. 

Cependant, les solutions de champ moyen ainsi obtenue ne sont pas uniques. Le choix 
d'un bon champ moyen repose essentiellement sur les proprieties du systeme considere, 
en particulier sur ses symetries. Ceci genere des difficultes majeures. Une quantite con- 
siderable de travail a ete fait sur ce point et il y a une litterature importante sur le 
sujet. En particulier, les systemes qui sont bien decrits par des modeles de type Heisen- 
berg sans frustration semblent , selon leur dimension, bien compris par l'introduction de 
l'etat ferromagnetique ou antiferromagnetique de Neel a temperature nulle T = E] ■ 
Cependant il n'en n'est pas de meme pour beaucoup de systemes de basse dimension- 
nalite (D < 2) et/ou frustres jHSEElllZl- Ces systemes presentent des caracteristiques 
specifiques. Une analyse extensive et une discussion, pour des dimensions D = 2, ont 
ete presentees par Wen [80J. La competition entre etat antiferromagnetique et liquide de 
spins ont fait l'objet de recentes investigations dans le cadre de la theorie quantique des 
champs a temperature nulle [7H1 Ell • 

La raison de ces comportements specifiques des systemes de basse dimension peut 
etre qualitativement relie au fait que les fluctuation quantiques et thermiques sont tres 
fortes, et detruisent ainsi l'ordre antiferromagnetique. Ceci motive une transcription de 
l'Hamiltonian en termes d'operateurs composes que nous appelons "diffuson" et "cooperon" 
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Chapter 1. Resume 



Le chapitre El compare le traitement exact (obtenue par la PFP) et moyen (par un 
multiplicateurs de Lagrange) de la contrainte d'occupation pour differentes approches de 
l'Hamiltonien de Heisenberg. 

En resume, nous montrons au chapitre El que la contrainte stricte d'occupation induit 
une difference quantitative de la temperature critique en la comparant avec les resultats 
obtenus par une contrainte d'occupation moyenne. En consequence il apparait un effet 
mesurable sur le comportement des parametres d'ordre. Cest la encore un point original 
de cette these qui a conduit a Particle donne en reference 

Avec l'occupation exacte de spin, la temperature de transition des etats de liquides de 
spins et antiferromagnetiques de Neel sont le double de la temperature critique obtenue 
par la methode des multiplicateurs de Lagrange. 

En revanche, la PFP ne peut pas etre employee dans une description en terme de 
"cooperons". Les cooperons sont a considerer comme des paires BCS et detruisent ainsi 
deux quasi-particules en faveur de la creation d'une nouvelle qui n'est autre qu'une paire 
du type BCS. II en resulte que dans ce cas le nombre de particules n'est pas conservee au 
contraire de la contrainte exacte. Or il se trouve que la PFP ne tolere aucune fluctuation 
du nombre de particule. 

Dans une description plus realiste nous devons tenir compte des contributions des 
fluctuations quantiques et thermiques qui peuvent etre d'une importance cruciale en par- 
ticulier au voisinage des points critiques. Le chapitre El se consacre a l'implication des 
fluctuations de phase autour du champ moyen dans l'etat 7r-flux. 

1.4 Du modele de Heisenberg a Paction QED3 pour 
des temperature finie 

Un ansatz de Neel n'est pas necessairement un bon candidat pour la description des 
systemes de spin quantique a deux dimensions. Nous avons montre au chapitre |U qu'un 
tel ansatz brise la symetrie SU(2) generant les bosons de Goldstone qui detruisent l'ordre 
de Neel. Un meilleur candidat semble etre le liquide de spins etant donne qu'il conserve la 
symetrie SU(2) intacte. Pour cette raison, le chapitre El se concentre sur la phase liquide 
de spins. 

L'Electrodynamique Quantique a deux dimensions d'espace et une de temps , la QED 3 , 
est un cadre commun qui peut etre utilise pour decrire les systemes fortement correles 
aussi bien que les phenomenes specifiques qui y sont relies comme la supraconductivite a 
haute temperature [23 12H EH1 EI] • Une formulation en theorie des champs du modele de 
Heisenberg en D = 2 dimensions d'espace fait correspondre faction initiale a une action 
QED 3 pour les spinons [2B1 EZ] • Avec cette description apparait le probleme du champ 
moyen et les questions correlees au confinement des charges tests qui peuvent conduire a 
l'impossibilite de determiner les contributions des fluctuations quantiques au travers de 
developpements en nombre de boucles EH1 E2] • 

Nous considerons dans le chapitre El l'etat 7r-flux initialement introduit par Affleck 
et Marston [21 Ell- L'occupation des spins par site du reseau est toujours fixee par 
la PFP. La PFP introduit un potentiel chimique imaginaire modifiant les frequences de 
Matsubara induisant en retour des modifications perceptibles au niveau de la temperature 
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de restauration de la symetrie chirale. La symetrie chirale est initialement brisee par la 
generation dynamique de masse des spinons. 

Nous nous concentrons sur la generation et le comportement de la masse des spinons 
due a la presence du champ de jauge U{1). Appelquist et cie. Ej ont montre qu'a 
temperature zero les fermions initialement sans masse peuvent acquerir une masse generee 
dynamiquement lorsque le nombre de saveur N est inferieur a la valeur critique N c = 
32/7T 2 . Plus tard Maris [52] a confirme l'existence d'une valeur critique N c ~ 3.3 au- 
dessous de laquelle la masse dynamique peut-etre generee. Etant donne que nous con- 
siderons des spins 5 = 1/2 nous avons iV = 2 et ainsi N < iV c . 

A temperature finie, Dorey et Mavromatos [21] et Lee [IB] ont montre que la masse 
generee dynamiquement s'annule pour une temperature T plus grande que la valeur cri- 
tique T c . 

Nous montrons que l'utilisation du potential chimique imaginaire introduit par la PFP 
[HE] modifie notablement le potentiel effectif entre deux particules chargees et double la 
temperature critique T c en accord avec les resultats presentes precedement [22J. 

Le potentiel chimique imaginaire reduit le phenomene d'ecrantage du potentiel d'inter- 
action entre des fermions test losrqu'on le compare a celui obtenu par la methode des 
multiplicateurs de Lagrange. 

Nous montrons egalement que la temperature de transition de la restauration de la 
symetrie "chirale", oil la masse des spinons m{(3) s'annule, est doublee par le potentiel 
chimique imaginaire de Popov-Fedotov. Le parametre d'ordre r = donne par Dorey 
et Mavromatos [21] et Lee [1H] s'en trouve reduit de moitie. La theorie, dans l'etat dans 
lequelle elle est actuellement, ne peut pas rendre compte de la mesure experimental de 
r qui est environ de 8 pour YBaCuO. 

Tous ces resultats ont donne lieu a Particle cite en reference [23]- 

1.5 Perspectives 

Marston [HH] a montre que pour retirer les configurations interdites de la jauge U(l) 
dans le modele antiferromagnetique de Heisenberg, un terme de Chern-Simons apparait 
naturellement et doit etre inclus dans Taction QED 3 . Cette contrainte supplementaire 
conduit a fixer a ir (modulo 2tt) le flux au travers de la plaquette formee par la maille 
elementaire du reseau de spins. De nouveaux travaux, que nous n'avons pas abordes 
dans ce manuscrit, montrent que la temperature de transition chirale ainsi que le rapport 
r = m/T c peuvent etre controles par le coefficient de Chern-Simons [20 . 

D'autre points interessants consernent la compactification du champ de jauge U(l) 
que nous avons utilise pour obtenir la masse dynamique des spinons. Dans le cas d'une 
theorie de jauge compacte des instantons apparaissent et interagissent avec la matiere 
(ici les spinons) et peuvent changer le comportement du systeme de spin jOH IBH] . Le 
confinement des spinons par les instantons reste un probleme ouvert. 
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Chapter 2 
Introduction 



After the discovery of high temperature superconductivity by Bednorz and Miiller JU] 
great theoretical efforts were devoted to find an explanation of the mechanism underlying 
the very high critical temperatures phenomenon. Many new superconductors were dis- 
covered amoung which one have cuprates with La 2 _ x Sr x CuOi and Bi 2 Sr 2 CaCu 2 0% are 
examples of hole doped superconductors, ruthenates Sr 2 Ru04, heavy-fermion metals such 
as UPts and ZrZn 2 , and organic materials like the well known k — (BEDT — TTF) 2 — 
Cu[N(CN) 2 ]Br. 

In particular cuprates behave differently from the conventional BCS superconductors. 
Experimentalists observed d-symmetry of the order parameter, strong electronic correla- 
tions and non-conventional but "universal" phase diagrams as shown in figure ETT1 [TTj . 



T 




Un- Under Optimaly Over 

doped doped doped doped 

Figure 2.1: Schematic phase diagram of a cuprate as a function of hole doping x and 
temperature T 

All cuprates share the same kind of atomic structure which consists of a layered struc- 
ture made of Cu0 2 -\&yeis which are considered as responsible for Cooper-pairing and 
sup er conduct i vity. 
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Chapter 2. Introduction 



In the following we concentrate on the C11O2 layers and more precisely on the anti- 
ferromagnetic insulating phase of the high-T c superconductors which corresponds to the 
undoped regime as shown in figure 12.11 The insulating phase of the cuprates like com- 
pound, the so called parent compound, can be modelled by Heisenberg models. 

The study of the insulating phase is motivated by the fact that under doping the 
parent compound should keep a memory of the correlations underlying the antiferromag- 
net insulating phase. In other words one believes that the superconducting phase must 
partially conserve some characters of the antiferromagnet insulating parent compound. 

Our work is devoted to the study of the effects of a strict site-occupation constraint 
caracterising the insulating phase in the description of quantum spin systems at finite 
temperature. This constraint consists in the enforcement of the occupation of each lattice 
site by exactly one S = 1/2 spin. The implementation of such a constraint was introduced 
by Popov and Fedotov 66J. Contrary to other methods which are based on the use of a 
Lagrange mutliplier [7J |H1 E] leading to an average occupations the Popov and Fedotov 
procedure (PFP) avoids this approximate treatment by means of the introduction of an 
imaginary chemical potential as we will show in chapter El We are interested to work 
out and analyse the application of the PFP onto antiferromagnet spin systems at finite 
temperature and to compare these results with the Lagrange multiplier method. 

Chapter El introduces the mathematical tools found in the literature and which will 
be used throughout this manuscript. We present the fermionization of spin Hamiltonian 
models and the PFP. We show that the PFP eliminates the unphysical Fock states in the 
fermionization of spin models. Finally we construct a path integral formulation of the 
partition function in imaginary time and show that the PFP induces a modification of 
the Matsubara frequencies which characterizes the fermionic propagator. 

In chapter SI we concentrate on a Neel mean-field phase description of quantum 
spin systems. The possible competition or phase transition between spin liquid states 
and an antiferromagnetic Neel state which may be expected to describe Heisenberg type 
systems in two dimensions where discussed recently [2E1 EH EH US] ■ We aim to study the 
pertinence of the Neel state ansatz as a mean-field approximation for finite temperature 
using the PFP. Quantum and thermal fluctuations contributions are worked out at the 
one-loop level taking the Popov and Fedotov procedure into account. The outcomes are 
compared with the spin-wave theory and extended to the XXZ-model by working out the 
magnetization and the magnetic susceptibility for temperature below and up to critical 
point. We discuss the degree of realism of the mean-field Neel ansatz. Indeed Ghaemi and 
Senthil [2H] have recently shown, with the help of a specific model, that a second order 
phase transition from a Neel mean-field to an ASL (algebraic spin liquid) may be at work 
depending on the strength of interaction parameter which enter the Hamiltonian of the 
system. 

The originality of my work stays here on the fact that magnetic properties have never 
been worked out using a strict site-occupation constraint. Chapter El remedies to it and 
an article on this point can be found in |21j . 

In chapter E] we consider different mean-fields possible choices of quantum spin sys- 
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terns at finite temperature in which each lattice site is occupied by exactly one spin 
imposed by the PFP. We also construct the formalism in which the occupation constraint 
is imposed on the average by means of the Lagrange mutliplier method with the aim of 
confronting these two approaches in the framework of Heisenberg-type models. 

The description of strongly interacting quantum spin systems at finite temperature 
generally relies on a saddle point procedure which is a zeroth order approximation of the 
partition function. The so generated mean-field solution is aimed to provide a qualitatively 
realistic approximation of the exact solution. The specific behaviour of low-dimensional 
systems is characterized by the fact that low- dimensionality induces strong quantum and 
thermal fluctuations, hence disorder which destroys the antiferromagnet order. This mo- 
tivates a transcription of the Hamiltonian in terms of composite operators which we call 
"diffusons" and "cooperons" which are the essence of the well known Resonant Valence 
Bond (RVB) spin liquid states proposed by Anderson 

The original point concerns the confrontation of the magnetization obtained through 
the PFP with the result obtained by means of an average projection procedure in the 
framework of the mean-field approach characterized by a Neel state. The same confronta- 
tion is performed for the order parameter which characterizes the system when its Hamil- 
tonian is written in terms of so called Abrikosov fermions (also called pseudo-fermions or 
spinons) [22]- 

Chapter El is devoted to a more realistic analysis in which we take care of the contri- 
butions of quantum fluctuations which may be of overwhelming importance particularly 
in the vicinity of critical points. 

Quantum Electrodynamics QED^+i) has attracted considerable interest in the last 
decade [SI EH ESI EEl EZJ since it is a common framework which can be used to describe 
strongly correlated systems such as quantum spin systems in 1 time and 2 space dimen- 
sions, as well as related specific phenomena like high-T c superconductivity [23 EE E01 EI] ■ 
A gauge field formulation of antiferromagnetic Heisenberg models in d = 2 space dimen- 
sions maps the initial action onto a QED3 action for spinons [2H1 EI] • We consider the 
7r-flux state approach introduced by Affleck and Marston [21 Ell- The strict site-occupation 
is introduced by the imaginary chemical potential proposed by Popov and Fedotov [HE! 
for SU (2) spin symmetry which modifies the Matsubara frequencies as will be explained 
in chapter El 

We show that at zero temperature the "chiral" symmetry is broken by the generation 
of a dynamical mass [SI E] which vanishes at finite temperature T larger than the critical 
one T c [2HEB1- 

We show that the imaginary chemical potential introduced by Popov and Fedotov 
[Up] modifies noticeably the effective potential between two charged particles, doubles the 
dynamical mass transition temperature T c and reduces the screening effect of this static 
potential between test fermions. This last point led to two articles given in [271] 

Chapter E] summarizes and comments our results and suggests further developments 
aimed to lead to clues to remaining open questions [2JJ 1221 EH] ■ 
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We present here a set of technical tools, extracted from literature and mainly taken 
from Negele and Orland's book jBI] , which are necessary in order to construct the partition 
function in terms of a functional integral. The fermionization of a Heisenberg spin- 1/2 
model is presented. The spin operators are described in the fermionic Fock space by means 
of the introduction of anticommuting creation and annihilation operators. In order to 
eliminate the unphysical Fock states which correspond to the presence of or 2 particles 
on each site of the system, we introduce the so-called Popov- Fedotov procedure (PFP). 
Finally we construct a path integral formulation of the partition function in imaginary 
time and show that the PFP induces a modification of the Matsubara frequencies which 
characterize the fermionic propagator. 



3.1 Fermionization of spin-1/2 Heisenberg models 
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3.1 Fermionization of spin-1/2 Heisenberg models 

Quantum antiferromagnet spin-1/2 models are of great interest for theoretical studies. 
Indeed they have deep connections with high-temperature superconducting materials such 
as Nd 2 - x Ce x CuO^ and La 2 - x (Sr or Ba) x CuO^ for which the undoped region of the phase 
diagram shows an antiferromagnetic phase described by a two dimensional square lattice 
of spin-1/2 particles, as explained in ^TJ IHOj- For example in the parent compound 
La 2 CuO± the copper is surrounded by six oxygen atoms forming CuO tetrahedron layers 
between La layers. Copper atoms interact with each other by means of a superexchange 
mechanism forming an effective two dimensional square lattice of spin-1/2 as shown in 
figures 13.11 




(a) (b) 

Figure 3.1: (a) Cristallographic structure of the cuprate superconductor La 2 _ x Ba x CuO± 
and (b) a Cu0 2 layer showing the (super)exchange coupling between the copper atoms. 

In the same direction Dagotto and Rice [IB] show that the two-leg spin-1/2 ladder 
materials, the vanadyl pyrophosphate (VO) 2 P 2 7 and the cuprate like SrCu 2 3 , can be 
modelised as a one-dimensional spin-1/2 Heisenberg model. 

In order to study these antiferromagnet spin-1/2 systems we will introduce the 
Heisenberg Anti-Ferromagnet Model (HAFM) with an external magnetic field Bi 

where the sums runs over the lattice sites r*j and r}. The effective (super)exchange coupling 
Jij (see A.Auerbach [H] chapter 1), acts between spins at position r*j and fj on a D- 
dimensional hypercubic lattice. For antiferromagnet systems Jy is negative and positive 
in the ferromagnetic case. The second term of the HAFM is the coupling between the 
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spins and the magnetic field B at each lattice site fj. The Lande factor and the Bohr 
magneton coefficients are absorbed in the magnetic field B. 

The spins S are vector operators, their components obey a SU{2) Lie algebra [S x , S y ] = 
iS z . Here we adopt the convention h = 1. The S = 1/2 spin vector operators can be 
expressed using the Abrikosov fermionic creation and annihilation operators f\ and f ia2 



Si — 2 fiffi ^°"i°"2 fio-2 (3-2) 

where o~i, a 2 =1,1 and the a aia2 vector components are Pauli matrices 

-(:;)■ -c^). ™ 

Explicitly spin operators at position i read 

St = flh (3-4) 

S- = flk (3-5) 
1 

2 1 

The creation and annihilation operators J %& and j\ verify the anticommutation relations 



St = oiflh-flh) (3-6) 



{ fiai 5 fja-2 } ~~ ^o-i ,cr 2 (3-7) 

{funJj**} = (3.8) 
ULfU = (3.9) 
Applying (|3.4jl . ()3.5j) and onto the physical Fock states |1, > and |0, 1 > one gets 

S+|0,1> = |l,0>=|t> (3.10) 

S-|1,0> = |0,1>=||> (3.11) 

5?|1,0> = -1/2|T> (3.12) 

S?|0,1> = 1/2| |> (3.13) 

The insertion of ()3.2|) into ()3.1|) generates the Fermionized Heisenberg Antiferromagnet 
Model which is quartic in the fermion operators. Using this formulation it is possible to 
write the Hamiltonian in different forms. Indeed the spin interaction term Si.Sj can be 
expressed as 

Si-Sj = Si.Sj 

- \t>\t> i ^ 

- 2 » lj+ 4 2 

+ ^ (3-14) 



2 
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where hi = fj a fi a , *D and C are quadratic in the fermionic creation and annihilation 
operators 

Va = 4/,t + ftfn (3-15) 

Cij = filfjl - filfjl (3.16) 

Neel states, Resonant Valence Bond (RVB) states jl] as well as 7r-flux states can be 
introduced by means of (|3.14|) in order to define different mean-fields, as done by using 
the large-N method developed by Affleck pQ, Affleck and Marston [2j , Read and Sachdev 

[53 EH EE. 

From equation ()3.6|) we see that the physically acceptable Fock states on a lattice 
site i are those for which the expectation value of the number operator either fj+fij or 

f}\fi,l is equal to one an the other equal to zero (corresponding to < S z >= ±1/2). The 
occupation by one fermion per lattice site is fulfilled if 

One way to implement this constraint consists of the introduction of a projector onto the 
Fock space using a Lagrange multiplier A 

P i = J vx^iiZrflf™- 1 ) (3.18) 

where A plays the role of a chemical potential in the framework of the path integral 
formulation. However such a propagator cannot impose a rigorous constraint ()3.17|) . The 
actual value of A is fixed by a saddle point method, that is to say mean-field equations 
will fix the mean value of the Lagrange parameter. In order to implement the constraint 
(13.17)1 in a strict way we shall introduce the Popov-Fedotov procedure. 

3.2 The Popov-Fedotov procedure (PFP) 

The introduction of this procedure is a key point in the present work. In the following 
we aim to study the consequences of its use in the description of spin systems at finite 
temperature. In the present section we show how it is able to enforce the constraint (|3.17|) . 

The Fock space constructed with the fermionic operators /, /' is not in bijective cor- 
respondence with the Hilbert space of spin states as shown in figure 13.21 Indeed, in Fock 
space and for spin-1/2 particles, the occupation of each site i can be characterized by the 
states Inj^n^ > with n i)(T 6 {0, 1} the eigenvalue of the occupation operator fj a fi a , that 
is the states |0, >, |1,0 >, |0,1 > and |1, 1 >. Since the sites have to be occupied by a 
single particle the unphysical states |0, > and |1, 1 > have to be eliminated. This is done 
by means of a projection procedure proposed by Popov and Fedotov [66J and generalized 
to SU(N) symmetry by Kiselev et al. |42j . 
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Hilbert space 

I T> 
I !> 



Fock space {|n-f,nj >} 
|1,0> 
|0, 1 > 
|0, > 

11,1 > 



Figure 3.2: Correspondance between the Hilbert space and the S = 1/2 Fock space 

We introduce the projection operator P = -^e l ^ N , where N = ^2fj a fi a is the number 
operator, into the partition function Z which then reads 



Z = Tr 



-0Hp 



(3.19) 



where H is the fermionized Hamiltonian of the systems and (3 = 1/T is the inverse 
temperature, with the convention ks = 1 for the Boltzmann constant. Define the trace 
of the operator O in Fock space as 



Tr[0]= [®< n i,hnu\ \0 (g)K T ,n a > (3.20) 

{m, a e [0,1]} \ i J \ i J 

where (S)| n i,T> n i,l > is the tensor product of the Fock states \riij, n^j, > for different lattice 

i 

site i. The contributions of the Hamiltonian should be equal to zero when applied to the 
unphysical states. Since Sf\unphysical >= S^\unphysical >= S-\unphysical >= all 
Hamiltonians build with only spin operators {Sf, S~ , S-} give an energy equal to zero 
when applied to unphysical states (H\unphysical >= 0). The action of Pj on each site j 
is such that the contributions of states 1 , >j and |1, 1 >j to the partition function Z 
eliminate each other. Indeed 

< 0, 0|, e-"* e*5*°|o, > j + < 1, 1|, e^ H J^\l, 1 
+ < l,0\ j e- f>H .e i *\l,Q> j + < 0,11^-^.6^10,1 >j 
= i (< 1,01^-^11,0 >j + < 0,l|j-e _/JH |0,l >j) 
Hence the partition function reads 



Z 



1 

7v 



.Tr 



(3.21) 



where Af is the total number of spin sites in the considered lattice. Equation (|3.21j) with 
the imaginary chemical potential [l = irk describes a system with strictly one fermion 
(spin- j or spin- j) per lattice site, in contrast with the usual method which introduces an 
average projection by means of a real Lagrange multiplier [HllZj. 

In the following the Heisenberg model will be studied in the path integral formulation 
using the Popov-Fedotov procedure. The consequences of its application will be compared 
to those obtained by means of Lagrange formulations. 
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3.3 Path integral formulation of the partition func- 
tion 

In the present section we construct the path integral formulation of the partition function 
(|3.21|) . First we shall define the coherent states of the fermionic Fock space. Then 
we derive from these states the appropriate properties which lead from the trace of an 
operator over Fock state into the trace over coherent states of the same Fock space. This 
leads to the functional integral expression of the partition function ()3.21j) . 

3.3.1 Fermionic coherent states 

The coherent states are an overcomplete linear combination of the set of states in Fock 
space. They are eigenstates of the annihilation operator, bosonic as well as fermionic [61 . 
Define |£ > as the coherent state of the fermionic Fock space 

oo 1 

l£>=E E (-l)-&-Gr-& , \m,-,n*,-,"k> (3-22) 

k=0 {n a }=0 

where \ni, . . . , n a , . . . , n k > is a Fock state constructed with the fermionic creation opera- 
tor yfij ■ ■ ■ {fl) na ■ ■ ■ (fk) applied on the vacuum state |0 >. The are Grassmann 
variables verifying the anticommutation relations 



= + = (3.23) 

= (3.24) 

£ = (£)* = (3.25) 

More properties on the Grassmann algebra are given in appendix El • 
Applying the annihilation operator f a on the coherent state |£ > we obtain 



oo 1 

f°\Z> = E E Ci 1 ---C a ---C k k fa\n l ,...,n a ,...,n k > 

k=0 {n a }=0 

_ i OL — 1 

OO 1 yy 

= E E er---e---4 n M-i)^ ni K---,n a -i,...,n,> 

fc=0 {^a} = 0, 
n a = l 

OO 1 "y 1 

= E E (-^^■■■C a+1 --^k k \^,---,n a ,...,n k > 

k=0 {n a }=0 

= UZ> (3-26) 

Thus ^ Q and |^ > are the eigenvalue and the eigenvector of the fermionic annihilation 
operator f a . Since fl increases the number of particles in any Fock state by one, |£ > 
cannot be an eigenvector of f\. 
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Coherent states are equivalently constructed combining the properties of the Grass- 
mann variables in exponential series such as 

|£> = e -S-^-|0>=IJ(l-e«/J)|0> (3-27) 

Q 

where |0 > is the vacuum state of the fermionic Fock space. Similary the bras of the 
coherent states |£ > verify 

<£| = < 0|e~ E « /a?s =< 0|e E «^ /a (3.28) 
<el/J = <£IC (3-29) 

Useful properties can be extracted from these basic equalities. 
3.3.2 Properties of the coherent states 

In this subsection we review some properties of the coherent states |£ > which will be 
used later in the construction of path integrals. The first one is the closure relation in the 
fermionic Fock space and the second is the trace of a fermionized operator O. 

We recall some main points concerning the closure relation of coherent states which 
reads 



J a 

A detailed demonstration can be found in the book of Negele and Orland [HI]. To proceed, 
we define the operator A as being the left hand side of equation ([3.30)1 . Following Negele 
and Orland, the first step in the demonstration of the equality ()3.30|) is to show that 
the operator A is proportional to the unit operator 1, the second step to show that the 
proportionality factor k is equal to one. 

Using the expression ()3.27j) and the definition of the derivation operator given in appendix 
|X]the commutation relation of fp and the operator |£ >< £| is given by 



[//?,!£ ><£|] = m >< £i - 1£ >< Z\fp 

= «><ei-ie>Jj<ei 



_d_ 



K><ei 



(3.31) 



Using the properties of the integration operator (see appendix EJ the commutation rela- 
tion of the operator A and the annihilation operator fp is equal to zero by virtue of the 
Grassmann algebra 
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<n* r \[fp,A]\nx> = < 



n-. 



n 



\i><i\\n x > 



«={A7.A} 

= (3.32) 

This equality can also be proven with the use of the creation operator /|. Therefore the 
operator A commutes with any operator composed of operators f a and /t. Since the 
Schur lemma stipulates that if an operator commutes with any operator then it must 
be proportional to the unit operator. A must be proportional to the unit operator I, 
A = k1. 

The matrix element of A between two vacuum states is given by 

<0|„4|0> = < 0| / Y[dC a d^ a e-^-^ a \^ >< f||0 > 

J a 

Y[de a du- Ea ^ 

a 

= JJl =< 0|kI|0 >= k (3.33) 

a 

Hence k = 1. 

The closure relation is very useful in order to define the trace of a fermionic operator 
O. Using the previously defined expression of the trace in the fermion Fock space given 
in Negele and Orland's book [HI] 



Tr [0\ = [Q$< n i,hnu\)0[(^)\n iih n i , l > 
= < n\0\n > 

n 

and inserting relation ()3.3()j) into the trace of the operator O, we obtain : 



Tr [O] 



Tr [O] 



a n 

a n 



(3.34) 



(3.35) 
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where < ({/„}, {/t}) |£ >= e -E.«e-0 ({-&}, {CD- 

These results complete the tools which enable the construction of the path integral 
(HOD). 

3.3.3 Partition function of many-body systems 

In this subsection we construct the partition function of a many-body Hamiltonian in a 
path integral formulation and we follow the prescription given in [HT]. We start with the 
grand-canonical partition function of a general fermionic Hamiltonian H in Fock space 



Z = Tr 



-I3(H-HN) 



(3.36) 



Using the Lie- Trotter relation lim (e uA -/ M e B / M ^ M = e equation ()3.36|) can be 



reexpressed as Z 



lim Tr 

M— »oo 



-s[H-nN 



M 



, with e = (3/M. Inserting the closure 



relation (|3.3U|) between each operator e £ ^ and using the expression of the trace 
over coherent states equation ()3.35j) . the partition function takes the form 



M 



lim / ]"] J! dC k ,J^a < -Ui\e-< H ^ N > >< £ 

J k=l a 

...ie 2 ><6ie- £ (^)|ei> 

. M 

^ i — 1 



M-l 



k=l a 



M 

£ 

fc=2 



+ e 



* r£a,l + £a,M 



+/^a,M} + #(c,i,-e 



(3.37) 



where a refers to the particle position and its spin, and k to the slice in which the term 
e -e\H-nN) a pp ear j n i\yQ Lie- Trotter formula. It is convenient to introduce the continuum 
notation £ a (r) to represent the set {£ a ,i, . . . , £ a ,k, ■ ■ ■ , £a,Af} and in the limit of M — > oo 
to define 



?«,* = ( r ) ^rs« (r) 



(3.38) 



The partition function Z of the many-body Hamiltonian H with a chemical potential /i 
takes the functional integral form 
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Z= f p^ e - ^ ^iE e(T)(3r-( 1 )( t ,(r)+//(e(r),( t ,(T))} 

■4»G8)=-€a(o) 

(3.39) 

M _ 

where = lim J] ]Q d£| a d£k a- Expression (j3.39|) is applicable for general fermionic 

M-*oo k=l 

Hamiltonian operators, in particular for the description of spin systems in terms of path 
integrals. In the next subsection we introduce a Fourier transform with respect to the 
imaginary time r. This leads to Matsubara frequencies which are shifted by the presence 
of the imaginary chemical potential \x. 

3.3.4 Modified Matsubara frequencies 

Here we show how the Matsubara frequencies are modified by the introduction of an 
imaginary chemical potential in the partition function of an spin-1/2 Heisenberg model. 
As we already saw above Popov and Fedotov introduced the imaginary chemical potential 
\x = ijp in order to remove the unphysical states of the Fock space as explained in section 

13.21 and in [HE]- Introducing this chemical potential in the term J(f^{EaG( r )(^ — 
A i )^a( r ) °f equation ()3.39|) the Fourier transform of £ Q (r) reads 



Ur) = Y.U^e*"* (3.40) 

Up 

and its reverse 



U"f) = \ f drU^e-^ (3.41) 
P Jo 

Here up = jf(n + 1/2) with n G 7L are the well known fermionic Matsubara frequencies. 
This leads to 



["drJ2C(r)(d T -^Ur) =DEGN^ " (3.42) 

Jo ~ „ Z P 



cop a 



Since the chemical potential is imaginary we see that we can redefine the fermionic Mat- 
subara frequencies by the change of variable 

tfF = "F-^ = j(n + l/4) (3.43) 
Redefining the Fourier transform of £ q (t) as 
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and shifting the partial derivative over r, ^ — [i — > the partition function takes the 
form 



The partition function itself does not change much by the introduction of the imagi- 
nary chemical potential /j, — if*. However the antiperiodic integration condition £ 7 (/3) = 
— £ 7 (0) needs to be changed into £ 7 (/3) = z£ 7 (0) by modification of the Matsubara fre- 
quencies. In the following chapters, we will show how these modified fermionic Matsubara 
frequencies up = jf(n + 1/4) change the behaviour of the physical properties of spin-1/2 
systems. 




(3.44) 




E C(r)w:Ur) + H({£(t)}, {&(t)}) (3.45) 
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Chapter 4. Mean-field and fluctuation contributions to the magnetic 
24 properties of Heisenberg models 



This chapter intends to present and discuss applications of the mean-field and loop 
expansion to the determination of physical properties of antiferromagnetic Heisenberg- 
type systems in spatial dimension D using the PFP. This original point of the present 
work was published in 1TT1 . 

Recent work on quantum spin systems discusses the possible existence of spin liquid 
states and in two space dimensions the competition or phase transition between spin 
liquid states and an antiferromagnetic Neel state which is naturally expected to describe 
Heisenberg type systems |2H1 EH HH HE] • It is also known that undoped superconducting 
systems show an antiferromagnetic phase [50J. 

In the following we focus our attention on a mean-field Neel phase description of 
quantum spin systems described by Heisenberg models. More precisely we present below a 
detailed study of the magnetization and the parallel magnetic susceptibility of Heisenberg 
antiferromagnetic spin- 1/2 systems on D-dimensional lattices at finite temperature. The 
aim of the work is the study of the physical pertinence of the Neel state ansatz, using the 
PFP, as a mean-field approximation in the temperature interval < T < T c where T c is 
the critical temperature [21]. In order to get a precise answer to this point we work out 
the quantum and thermal fluctuation contributions beyond the mean-field approximation 
under the constraint of strict single site-occupancy EU |22l E21 EE] which allows to avoid 
a Lagrange multiplier approximation [Jj. The results are also extended to anisotropic 
XX Z systems and compared to those obtained in the framework of the spin-wave theory. 
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4.1 Nearest-Neighbour Heisenberg model 

In order to focus on the essential properties of a Neel state mean-field we consider a 
Heisenberg model which present a bipartite spin lattice. A bipartite lattice can be split 
into two disjoint sublattices A and B, where J^- connects only i G A to j G B as defined 
in [H]. The simplest isotropic Heisenberg model showing a bipartite lattice reads 



H = --^JijSi.Sj + ^Bi.Si 



(4.1) 



with 



ffe{Si,...,a D } 



(4.2) 



where J is the negative antiferromagnet exchange coupling working between nearest neigh- 
bour sites % and j separated by the lattice vector ff e {d\, . . . , do} on a D-dimensional 
lattice. Keeping the static external magnetic field always fixed in the Oz direction the 
Hamiltonian reads 



(4.3) 



<i,3> 



where the sum ^ runs over the nearest neighbour site < i,j > at position r. \ G A and 

<i,j> 

fj G B. 



Figure 4.1: A two dimensional bipartite lattice system in a Neel state. 



Figure l4~D is a two-dimensional representation of a bipartite lattice. The spin sublattice 
A and B are represented by blue and red arrows, dotted lines materialize the exchange 
interaction J < between spins (the blue and red arrows). 



4.2 Spin- wave theory 



A few years ago Coldea et al. [Tj\ measured the magnetic excitations of the square-lattice 
spin-1/2 antiferromgnet and high-T c parent compound La^CuO^. They showed that the 
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inclusion of some interactions beyond the nearest-neighbour Heisenberg term (|4.2JI leads to 
a good description of the dispersion relation observed by high-resolution inelastic neutron 
scattering in the framework of spin-wave theory. 

A modified spin-wave theory developed by Takahashi [77j led to the Auerbach and 
Arovas equations 7\ which were obtained by Schwinger-boson formulation. 

In this section we calculate the sublattice magnetization as well as the magnetic sus- 
ceptibility in the framework of spin-wave theory applied on the nearest neighbour antifer- 
romagnet Heisenberg model (|4.3j) . Later these results will be compared to those obtained 
using the Popov- Fedotov procedure with the Hamiltonian (|4.3|) . 

4.2.1 The Holstein-Primakoff approach 

In a broken symmetry phase such as the Neel state at least one spin component shows a 
non-zero expectation value. For small temperatures the fluctuations about this expecta- 
tion value can be studied by means of the Holstein and Primakoff (H-P) spin-deviation 
creation and annihilation boson operators dj and a\. 

Since the external magnetic field is applied in the Oz direction the non-zero spin 
component is S z . Following H-P we express the spin operators in the form (HE] 



S+ = ^S-nfOi 

Sf = S-nt (4.4) 

where = a\di = S — Sf is the so called spin-deviation operator, eij and a\ are boson 
operators. Here 5 is defined from the relation 5 2 = 5 (5 + 1) of the quantum spin vector 
5. For small temperatures and/or for large-5 the expectation value of the spin-deviation 
operator < n? > is small compared to 25 and leads to 



+ (4 ' 5) 

Since the Hamiltonian (|4.3|) shows a bipartite structure on the lattice the spin expectation 
values are of opposite sign from sublattice A to sublattice B as depicted in figure 14.11 
We can define two type of spin operators depending to which sublattice (A or B) they 
belong. Introducing (j4.5)) into ()4.4j) and neglecting all terms nf/S the spin components 
for sublattice A read 



Q+ 

D A,i 


= v25aj 




= a\V2S 


QZ 


= S-nt 



(4.6) 



and for the sublattice B 



4.2 Spin- wave theory 



27 





= V2Sa] 


S B,j 


= a,j V2S 


qZ 


= S + n« 



Following Igarashi [SB], Kubo [IB], Oguchi JHSJ and Takahashi [ZZ] the use of (|4.fij) and 
fj4~7j) in (jOJ) leads to 



H. 



sw 



ieA 



fj — fj + rf 

+ ^2 B i S A,i + B i S B,j 

ieA jeB 



-^z\J\S 2 + z\J\S 



ieA 



jeB 



+\j\s^2^2 { a i a j + a ! a i) 



ieA jeB 



(4- 



where M is the total number of spin-S 1 in the system and z = 2D is the coordination 
of a spin in a D-dimensional hypercubical lattice. The magnetic field Bi is redefined as 
Bi = d§fg- Since we admit that the expectation value of the spin-derivation operator 
< n a > is small only quadratic terms in the boson operator a appear in the Hamiltonian 
()4.8|) which amounts to neglect the interaction between spin-waves. 
The spin-wave partition function Zsw a t finite temperature reads 



Z sw = Tr [e^Hswm})] (4.9) 

Magnetization and susceptibility of the spin system will be extracted from ()4.9|) taking 
derivatives with respect to the magnetic field Bi. 



4.2.2 Neel Spin-wave magnetization 

In order to derive the sublattice magnetization rriA(B) from the spin- wave free energy 
Fsw — —jj^Zsw we set Bi = B if i belongs to sublattice A sites and Bi = —B if it 
belongs to sublattice B. Taking the derivative with respect to B in (|4.9jl we obtain 



m A 



J7\l2 < ^Xi >SW < S B,j >SW j 

\ieA jeB J 



1 d \ 1 d 

NpoB \b=o N oB b=o 



sw 



(4.10) 
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Define the Fourier transform of the boson operator a 



-7^= E <> 



(1) ik.fi 
k C 



keSBZ 



^= E » 



(2) ik.fj 
k C 



(4.11) 
(4.12) 



keSBZ 



where SBZ is the Spin Brillouin Zone which is shown in Figure l4~2l (shaded area) inside 
the lattice Brillouin Zone of a two dimensional bipartite spin system (large square). More 
details on the construction of the Spin Brillouin Zone are given in appendix iBl 

k,. 




Figure 4.2: Two dimensional Spin Brillouin Zone (shaded area) 



Fourier transforming (|4.8J) with the definitions (|4.11|) and (|4.12j) the Hamiltonian goes 
over to 



H, 



sw 



\J\S 2 +Afz\J\BS 2 + z\J\S y \(l-B)bfhf + {l-B)bfh 



+z\J\S J2 7^ 



keSBZ 
k — k' k — k 



(4.13) 



keSBZ 



where 7^ is the Fourier transform of XE<5(r*j ~ + V) anc ^ V was defined in equation 
(14.2)1 . section |4~T1 7 is related to ff by 



Ik 



z 



E 



Jk.if 



(4.14) 



??e{±ai,...,±a D } 



Notice that the coordination z is simply related to the lattice vectors by z = Yl 1. 

ffe{±ai,...,±a D } 

The Hamiltonian (J4.13J1 can be diagonalized by means of a Bogoliubov transformation 
applied to the boson operators b, [HJ [T3] 



6« 



(+),* 



(4.15) 



(-),-* 
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where and t>g are real coefficients verifying ui — vi = 1 as a consequence of the 
commutation relations obeyed by the boson operators g and {3,-. g. We set = 
cosh and v? = sinh 9?. Using 1)4.15)1 leads to the diagonalized Hamiltonian (J4.1H)) 



H sw = -^-z\J\S 2 + z\J\AfBS 2 -^-z\J\S(l-B) 



with 



tanh2% = -j^B) (4-17) 



z^S.Ul-B)-^ (4.18) 



Here cj^ is the spin-wave spectrum and £^(±) fc * s ^ ne num b er operator of magnons 
occupying the energy "level" Ut- In the absence of the magnetic field B and for k close 
to Q = or 7? (7? = (7T, ...,7r)) the dispersion relation shows a relativistic spectrum 
uj: ^ c\k — Q\ where c is the spin- wave velocity for a D-dimensional hypercubical lattice, 

c=vb\j\s. 

The sublattice magnetization rriA = —itlb is derived from the spin- wave free energy 
J= sw = -i In Tr [e~^ Hs ^] and reads 



1 



m A 



N tanhfwr' f\ ~2 
keSBZ 2 k V ^ 



keBZ 



1 



^V" J+ 2JVT^ <419) 



— — is the boson occupation number. Notice that in the second equation of rriA 

K e fe — 1 

the sum runs over the whole Brillouin Zone. Later the magnetization ()4.19|) will be used 
as a reference for comparaison with the magnetization obtained by means of our method 
using the PFP. 



4.2.3 Neel Spin-wave susceptibility 

With a uniform magnetic field Bi = B applied on the spin system the spin-wave suscep- 
tibility is obtained from the free energy by 



XII 



d 2 T, 
OB 2 



sw 



B=0 



(4.20) 



where J-< 



sw 



- \ In Tr [e ^ SHF ] and the spin- wave Hamiltonian Hsw 
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H 



sw 



z\J\S 2 + z\J\S 



£ (1 - B) a}* + £ (1 + B) aja. 



jeB 



+ (a^j + a] a] 



(4.21) 



i&Aj&B 

A Bogoliubov transformation leads to the diagonalized spin-wave Hamiltonian 



H SW = ~jz\J\S 2 + z\J\S + 

keSBZ 



+ H - b) ^A-A ( 4 - 22 ) 

tanh(2%) = - Tk (4.23) 



z\J\S^l^l (4.24) 
B = B/z\J\S. The free energy is then given by 



sw 



-^liaZsw 



-y 

R ^ 







Put - In 



1 -e 



1 -e 



(4.25) 



keSBZ 

Finally the spin-wave susceptibility reads 



d 2 Fsw I 
X{1 8B 2 \b=o 

= 2(3 J2 n A n k + 1 ) (4-26) 

keSBZ 

where = — is the boson occupation number of the magnons. 



4.3 Effective action 

After this review on spin-wave theory we return to the study of the Heisenberg model by 
functional integrals. The aim of this work is to extract the magnetization and the spin 
susceptibility from the path integral formulation (see section equation (|3.45jl ) for the 
Hamiltonian f)4.3|) of a bipartite spin system in a Neel state. As we shall see later the 
analytic development of the effective action necessitates some approximations. In order to 
be able to appreciate the pertinence of these approximations, that is the Neel mean-field 
with the one-loop corrections for example, we define a "Ginzburg-Landau parameter" 
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which evaluates the relative importance of the quantum and thermal fluctuations on the 
mean-field Neel state. We compare the magnetization and the susceptibility worked out 
from functional integral formulation with those obtained from the spin-wave theory for 
the Hamiltonian (|4.3j) . 

4.3.1 The Hubbard- Stratonovich transform 

As was shown in chapter |3] the Hamiltonian ()4.3|) can be expressed in terms of creation 
and annihilation fermion operators /' and / by replacing the spin operator Si by ([3.2)1 . 
Performing this transformation and injecting the resulting fermionic Hamiltonian ()4.3|) 
into the path integral ()3.45|) we see that the term Si.Sj in the exponent leads to a quartic 
expression in the Grassmann variables 

i 

- E J te^K^v^)) • (&^)*^ >4 (t)) (4.27) 

<i,j> 

where repeated indices mean summations over them. Recall that a are the 577(2) Pauli 
matrices defined in 13.21 Since integration over the Grassmann variables of the path 
integral (|3.45|) is not possible as it stands we have to reduce the action (|4.27J) to a quadratic 
expression by means of a Hubbard- Stratonovich (HS) transformation [HUES]- We define 
the Hubbard-Stratonovich auxiliary field action So as 



So 



if (Tj 



where (J -1 )^ is the inverse of the coupling matrix and is the Hubbard-Stratonovich 
auxiliary field. The matrix (J" 1 )^. always exists if one considers periodic bondary condi- 
tions on the spin system. Adding the HS action Sq to the fermionic Hamiltonian H (£*, £) 
and performing the change of variable 

^(r)=^(r) + ^4(r) (4.29) 



where Sj(t) = ^ >ai {T)a aic ^j ja2 {T) leads to the Hubbard-Stratonovich transformed Hamil- 
tonian which reads 



Up (Tj 



+ H (£*, - S [uj(t)} + (&{t) + 4.) Mr) (4.30) 
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and the partition function 



Z 



z n 



J <p(J3)=<p{p) ' •/;„.,( »')_,;„, (ih 

d 



(4.31) 



5[^r,C] = ^a(r)^(r) + 5 b(r)] + ^(^(r) + J B l ).S ? J (r) (4.32) 
Here Z stands for the partition function of the HS auxiliary field and reads 



V(pe~ f° drS ° mr)] 



(4.33) 



M 



and Dip = lim |~[ Wdtp* k dtp\ k difl k as explained for the Grassmann variables in subsec- 

M^oo- - ■ ' ' ' 

tion 13. 3. 31 



k=l i 



4.3.2 Integration over the Grassmann variables 

The action ()4.32|) can also be written as 

S[^C^} = SoMr)-B] + J2(^ ) M ^(|'l(r) ) (434) 

after the variable shift <Pi(r) — > 0i(r) — B, L . The matrix Mj(r) contains the factor of the 
quadratic terms in the Grassmann variables 



MAt) 



KO0 £ - ^f(r) 



(4.35) 



Fourier transforming (|4.34j) with the definitions of subsection 13.3.41 one gets 



0i{-UB,n) - Bi (fj(uJB,n) ~ Bj 

^B,n 1,3 

+ pE E UW 3 **) ) m(p - g) ( | T [~H 



(4.36) 



where Wb,»i = are boson Matsubara frequencies since from equation (|4.29|) the aux- 
iliary field ip needs to be periodic in r so that (pi(/3) = <Pi(0). The modified fermionic 
Matsubara frequencies 0p )P = %(p+l/4) are defined in subsection !3.3.4l After the Fourier 
transform the matrix M,- reads 
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Mi(p - q) 



(4.37) 



i^F,pSp, q + ±iPi(Cj F> p - UJ F ,q) {^F,p ~ &F,q) 

\<Pt{&F,p ~ & F , q ) i0JF,p5 p ,q ~ (&F,p ~ &F,q) 

Then after integration over the Grassmann variables using equation (jA.24j) one obtains 



2=1- f V<pe- S °ff [ ® 



(4.38) 



where the effective action S e ff reads 



rP 

S eff = / S [(p(r)-B}-y2\ndet/3M i 
Jo 



(4.39) 



4.4 Mean-field equation and One-loop contributions 



The mean-field equations are obtained from 



ss, 



Sip 



which is the stationnarity condi- 



tion in the application of the least action principle. Assuming that the mean-field and 
fluctuations contributions of the effective action S e ff can be identified and separated, the 
matrix M can be decomposed into a mean-field part Gq 1 and a fluctuation contribution 
Mi and thus reads 



M = -G l + M l (4.40) 

where the mean-field matrix Gq 1 depend on the choice of the mean-field Hubbard- 
Stratonovich auxiliary field (p and M\ is composed of the fluctuations 5ip = 0—0- At first 
glance it does not seem obvious to choose the mean-field components of the HS auxiliary 
field (p. From 1)4.29)1 it is clear that the mean- field of (p is related to the mean- field of the 
spins and thus <p should show the same symmetries as the spin mean-field. When the tem- 
perature is increased it is expected that thermal fluctuations Sip become more and more 
important and the auxiliary field (p may move away from the mean-field (p. Considering 
the temperature limit T — > the boson (and also fermion) Matsubara frequencies UB, n go 
to zero for any value of n. Then the relevant mean-field Fourier components of Lp{uoB, n ) 
are those for which ujs, n = 0. Fourier transforming M with respect to the imaginary time 
t and extracting the mean-field part (p from M one obtains 



G 



Op,q 



1 



det Gr 



[iu FtP - ±(p?(uj FjP - U Fjq )] 5 p>l 



UF,q)8p,( 



1 



det G„ 2 



\ip i (u F>p ~ U F ,q)5 p , q 



1 



det G r 



[iOJ FtP + ~(pf(LU FtP - LO F>q )] S PtQ 



(4.41) 



and the fluctuating part 5ip are given by 
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M 1 



p,q 



\8(p?(w FiP - u Fq ) \6(p { {uj Fp - u Fq ) 



with 50i(u Fp — u Fl 



;S(ff(uJ Fp - U) Fq ) -^5(fl(uJ F p - UJ Fq )_ 

0i(u Fp - u Fq ) - 0i(uJ Fp - uJ Fq )5 m and det G p 



(4.42) 



2 J 



. The expression lndet(/3M) in the effective action (|4.39|) 



can now be developed into a series 



lndet(/3M) = lndet/3 [-Go^f -G M 1; 

= In det (-/5Gq *) +Tr ln(l - G M X ) 



oo 1 

lndet(-/?G x ) - Tr{^ -(G Mi) n } 



71=1 



The first term lndet(— f3G a l ) leads to the expression ^ ln2 cosh ^\\0i(u>B = 0)|| and the 
effective action over the auxiliary field <p reads 



S, 



eff 



p/3 a T °° 1 

/ drS o [0(r)]- Vln2cosh^||^(^ = 0)||+Tr V-(G Mi) n 

1/0 t n=l 



(4.43) 



The first term n = 1 in the sum over n gives the contributions at the first order in the 
fluctuations 8<p, the second one the one-loop correction to the mean-field for n = 2. It is 
quadratic in Sip. Hence in a loop expansion beyond the mean-field approximation the 
effective action given by ()4.43|) is a Taylor series expansion in powers of Sip. To second 

order (one-loop contribution) in the fluctuations 8<p of = ip + 80. 



S, 



eff 



s. 



eff 



+ 



ss, 



eff 



[0] 80 [0] 



8<p 



18 2 S, 



eff 



2 80 2 



80 2 + O(8p 



(4.44) 



We now give a more precise definition of our mean-field (p. The partition function 

80 = 



(|4.38|) can be worked out by means of a saddle-point method in which g ^ 
and leads to the equation 



-s, 



Z = e 



eff 



V8ipe 



1 °eff 
2 5(32 



where the mean-field solutions verify the self-consistent set of equations 



(4.45) 



E 

3 



(■I % 0i-Bj 



2 pi 



tanh 



~P<Pi 



(4.46) 
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These equations lead directly to the cancellation of the first order term in Sip as worked 
out in appendix O 

In the following we consider a Neel mean- field order (pi(r) = (—l) 7T ' ri (f z e z = <p z e z 
where tt is the Brillouin spin sublattice vector as defined in section 14.21 A magnetic field 
aligned along the direction e z is applied to the system. The partition function can be 
decomposed into a product of three terms 



Z — Zmf-Zzz.Z^ 
where Zmf, Z zz and Z + _ are given by 



(4.47) 



Zmf 

Z zz 

Z+_ 



~ S eff 

e 



^0 



Z n 



1 gjgg/j 

2 a v z2 



V5<p z e 



5ip z 



V(6<p + ,5<p~)e 



2 9y>+ dip 



8ip~*~ .Sip 



(4.48) 
(4.49) 

(4.50) 



with 



S, 



eff 



M 2 z - 



(4.51) 



The one-loop corrections (^y^) 2 and <5y9 + 5<^ terms are worked out in details in appendix 
lUland read 



1 S 2 S, 



eff 



in, 



^tanh' = 



^(-w B )(5^(a; s ) 



1 d 2 S, 



2 <fy+<fy- 



=E£§ 



5 (n, 



l tanh(f^) 

2 iff- iuj B 

5tpj{uj B )5<pj{-uj B ) 



5(pf{-uj B )5(p j (u B ) 



(4.52) 



Zmf is the mean-field contribution, .2^ and are the one-loop contributions respec- 
tively for the longitudinal part 5ip z and the transverse parts of ip, 5tp + ~ , which take 
account of the fluctuations around the mean- field value (p z . 
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The contributions Z zz and are quadratic in the field variables 8(p z , 5ip + ~ and 
can be worked out in the presence of a staggered magnetic field Bf. Studies involving 
a uniform magnetic field acting on antiferromagnet quantum spin systems can also be 
found in ref. 42 . 



4.5 Magnetization and susceptibility for D-dimension- 
al systems 

4.5.1 Relation between the Hubbard- Stratonovich auxiliary field 
and the magnetization 

We define by J~mf = — -q]rZmf the mean- field free energy where Zmf is given by 
equation ()4.48j) . The local fields {<fii} can be related to the local magnetizations {fhi}. 
Using fhi = — 9 ^ F one gets 



MF 



1 d 



s, 



(4.53) 



and from the mean-field relation (j4.46|) one deduces also the relation 



in; = tanh — u)i 

2<pi 2^ 



(4.54) 



where tpi = \\ipi\\. Considering the Neel state tpi = <fie z and keeping the external magnetic 
field Bi applied in the direction Oz one gets 



ft - B t 



Z, _1 _ 

—tanh 2rrii 



rrij 



(4.55) 
(4.56) 



Combining (|4.55jl and (|4.56jl the self-consistent mean-field equation of the magnetization 
fhi is obtained by means of the set of equations 







tanh 1 2rrii = Bi + J l 



rrii 



1 

- tanh 
2 



(4.57) 
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4.5.2 Linear response theory 

The magnetization mean-field equation 1)4.57)1 can be solved in the linear response theory. 
If the applied magnetic field is weak enough the mean-field magnetization fhi can be 
developed linearly with respect to the magnetic field 



rhi = (-l) rV m + AmBi + 0(B 2 ) (4.58) 

where tt is the Brillouin vector coming from the existence of the sublattices in the Neel 
state defined below f)4.46)) . the lattice position, m the sublattice magnetization and 
Am the linear coefficient in the magnetic field Bi. According to the dependence of the 
magnetic field on the site i of the lattice and by inspection of ()4.58j) and ()4.57|) the linear 
coefficient Am reads 




£(l-4m 2 ) 



4 



, when Bi = (-lY^B 

i + |D|j| (1 -4.) ,whenB, = B. 

where m is the mean-field sublattice magnetization without external magnetic field which 
verifies 



m = -th-D\J\m (4.60) 

The solution depends on the dimension D and the coupling constant J of the nearest- 
neighbour Heisenberg model ()4.3|) . The self-consistent equation ()4.60|) is easily solved 
numerically using the Newton method [HZj- As can be seen on figure 1431 the mean-field 
magnetization saturates at 1/2 for the temperature T = and vanishes at the critical 
temperature T c = D\J\/2. Depending on the configuration of the magnetic field one can 
either compute the magnetization with Bi = (— ly^B or the susceptibility with Bi = B 
as was previoulsy explained in section 14.21 concerning spin- wave theory. 



4.5.3 Neel magnetization with fluctuation corrections : results 

Substituting equation 1)4.58)1 with the magnetic field configuration Bi = (—l) Ti - n B and 
relations (j435]) . flOg} into equations (fQ5|) . (jO^ and f)43U|) integrating over the HS 
auxiliary fluctuation field 5ip one obtains the free energy T . The derivation of the free 
energy is given in details in appendix O The components J-mf,J~zz an d F+- of the free 
energy for a linear approximation in the magnetic field read 



J~mf = ND\J\ (m + Am .B) 2 - ^ In cosh ( ^ [B + 2D\J\(m + A 



rho.B)} 



(4.61) 
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(4.62) 
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sinh ( § ([B + 2D\J\{m + Arh .B)} 2 - [2D\ J\^(fh + Am .B 



1/2 



sinh (f [B + 2D\J\(m + Am .B)]) 

(4.63) 

The magnetization m on site i is the sum of a mean- field contribution m = — \ 9ln ^ F , 
a transverse contribution <5m + _ = —■^ 91 q^~ an< ^ a longitudinal contribution 5m 22 = 
p. For a small magnetic field B a linear approximation leads to m = rh + <5m^ + 
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P dB 

5m + _ where 
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keSBZ 
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(4.64) 
(4.65) 



A/" is the number of spin-1/2 sites, Am 



; 1-4. m 2 



(4.66) 

and 7£ = -^J^ cos(k.ff) as 

2 — I- iv- I fj 

defined in section l4~2l 

At low temperature (T — > 0) it is seen that the magnetization goes over to the corre- 
sponding spin- wave expression [SHI EU EEJ E] and reads 



l-fD|J|(l-4.m 2 ) 



m 



t6sBztaIlh( «M x/1 _ 7l , n_ 7 | 



m + Sm 



(4.67) 



where m = 1/2 is the mean- field contribution and Sm is generated by thermal and quan- 
tum fluctuations. 
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Figure 4.3: Magnetization in a 3D Heisenberg antiferromagnet cubic lattice. Dotted line : 
Mean-field magnetization, Dots : Spin-wave magnetization, Full line : One-loop corrected 
magnetization. 

Figure 14.31 shows the magnetization m in the mean-field, the one-loop and the spin- 
wave approach (as given in section fO|) for temperatures T < T c where T c = D\J\/2 
corresponds to a critical point. One observes a sizable contribution of the quantum and 
thermal fluctuations generated by the loop contribution over the whole range of tempera- 
tures as well as an excellent and expected agreement between the quantum corrected and 
the spin-wave result at very low temperatures. 

The magnetization shows a singularity in the neighbourhood of the critical point. 
This behaviour can be read from the analytical expressions of <5m + _ and 8m zz and is 
generated by the |fc| = mode which leads to 7^ = 1 and by cancellation of m. The Neel 
state mean-field approximation is a realistic description at very low T. With increasing 
temperature this is no longer the case. The chosen ansatz breaks a symmetry whose effect 
is amplified as the temperature increases and leads to the well-known divergence disease 
observed close to T c . Hence if higher order contributions in the loop expansion cannot 
cure the singularity the Neel state antiferromagnetic ansatz does not describe the physical 
symmetries of the system at the mean-field level at temperatures in the neighbourhood 
of the critical point. Consequently it is not a pertinent mean-field approximation for the 
description of the system. 

The discrepancy can be quantified by means of the quantity where Am = m—m = 
5m zz + <5m + _. Figure l4~4l shows the result. The relation < 1 (Ginzburg criterion) 
fixes a limit temperature Tu m above which the quantum and thermal fluctuations generate 
larger contributions than the mean-field. For 3-D systems this leads to Tu m ~ 0.8T C , for 
2D systems the criterion is never satisfied except maybe for very low temperature , see 
figure 14.51 

The pathology is the stronger the smaller the space dimensionality. It is also easy 
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Figure 4.4: Ginzburg criterion =^ for the 3D Heisenberg model. 



to see on the expression of the magnetization that, as expected, the contributions of the 
quantum fluctuations decrease with increasing D. As can be seen in the figure 14. 5[ the 
saddle point breaks down earlier in two than in three dimensions. 

In fact, the Heisenberg model spin- wave spectrum shows a Goldstone mode as a con- 
sequence of the symmetry breaking by the Neel state. When \k\ goes to zero 



ZDSJI-^ 



lim LJr 

k^O 



(4.68) 
(4.69) 



The zero mode destroys the long range order in ID and 2D as expected from the Mermin- 
Wagner theorem 



4.5.4 The susceptibility : results 

We consider the parallel susceptibility x\\ which characterizes a magnetic system on which 
a uniform magnetic field is applied in the Oz direction. The expression of x\\ decomposes 
again into three contributions 



1 d 2 T 
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Xmf + Xzz + X+- 



(4.70) 
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Figure 4.5: Comparison of the Ginzburg criterion applied to a 2D (dotted line) and 3D 
(full line) Heisenberg model. 
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(4.73) 
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The behaviour of \\\ is shown in Figure l4~o1 where we compare the mean-field, spin-wave 
and the one-loop corrected contributions for a system on a 3 D cubic lattice. One observes 
again a good agreement between the quantum corrected and the spin wave expressions at 
low temperatures. For higher temperatures the curves depart from each other as expected. 
The mean-field contribution remains in qualitative agreement with the total contribution. 
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Figure 4.6: Parallel magnetic susceptibility at 3D for the Heisenberg model. Dots : 
Spin- wave susceptibility. Dotted line : Mean- field susceptibility X\\ MF - F uu li ne : total 
susceptibility (x\\ MF + *X||) - 



4.6 The XXZ-model 



To shed light on the fluctuations created by the symmetry breaking on an Heisenberg 
model when a Neel state is used as a mean- field state we add an anisotropic term — ^(1 + 
5)S*Sj. In this case we consider a so called XXZ-model. The corresponding Hamiltonian 
of the system can be written 



H xxz = -^J2 ( S ? S * + S i S J + + S ) S " S i) ( 4 - 76 ) 

<ij> 

where 6 governs the anisotropy. The self-consistent mean-field magnetization of the XXZ- 
model reads 



m = -tanh^D\J(l + 8)\m (4.77) 

Following the same steps as for the Heisenberg model we can compute the linear response 
of the system to the application of a magnetic field Bi in the direction of the Oz axis. 
The corresponding response to the mean-field magnetization m, = (— l) l fn + Arhxxz-Bi 
reads 



f (1 -4.m 2 ) 

Amxxz = s — — 4.78 

1 -fD|J(l + 5)| (l-4.m 2 ) V ' 

The magnetization is derived from the free energy given in the appendix and reads 
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m = m + 5m zz + <5m + _ 



(4.79) 



where 
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(4.82) 



tanh (Wlmy/l-^) 2 ) 



(4.83) 



The critical temperature Tf xz of the mean-field magnetization for the XXZ-model reads 



T. 



XX z 



D\J(1 + S)\ 



(4.84) 



For the zero temperature the excitation spectrum of the spin-wave obtained from the 
previous magnetization expressions leads now to a finite \k\ — energy 



u % = ZDS\ll 



1 + 5 



Ik 



(4.85) 
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By examination the expressions show that the zero momentum mode is no longer 
responsible for a breakdown of the saddle point procedure near 

T xxz = ^J(i±M However 

the magnetization of the XXZ-model remains infinite near Tf xz . This is due to the 
common disease shared with the Heisenberg model that the mean-field magnetization 
appearing in the denominator of 5m + _ goes to zero near the critical temperature. One 
sees that the mean-field Neel state solution makes only sense at low temperatures, that 
is for T < Tu m , whatever the degree of symmetry breaking induced by the mean-field 
ansatz. 

One notices that the spectrum no longer vanishes in the limit k — > when the 
anisotropic coupling 5^0. Since this is the case the transverse modes need a finite 
amount of energy to get excited. Thus they are not Goldstone modes by virtue of the 
Goldstone theorem [3122] • The Neel state does not break the 0(2) symmetry but only 
the discrete one Z2 (in the anisotropy direction Oz) of the XXZ-model. 
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4.7 Summary and conclusions 

In the present chapter we aimed to work out the expression of physical observables (mag- 
netization and susceptibility) starting from a specific mean-field ansatz and including 
contribution up to first order in a loop expansion in order to investigate the effect of fluc- 
tuation corrections to mean-field contributions at gaussian approximation. The mean-field 
was chosen as a Neel state which is an a priori reasonable choice for spin systems de- 
scribed in terms of unfrustrated bipartite Heisenberg model. The results were compared 
to those obtained in the framework of spin-wave theory. The PFP does not qualitatively 
change the well known magnetic properties of the Heisenberg model. However it must 
be pointed out that strict site-occupation had never been taken into account. This work 
remedies to this lack and was published in [21]. 

The number of particles per site was fixed by means of an strict constraint implemented 
in the partition function. It has been shown elsewhere [HI [22] that this fact introduces 
a large shift of the critical temperature compared to the case where the constraint is 
generated through an ordinary Lagrange multiplier term. 

At low temperature the magnetization and the magnetic susceptibility are close to the 
spin- wave value as expected, also in agreement with former work Quantum correc- 
tions are sizable at low temperatures. With increasing temperature increasing thermal 
fluctuations add up to the quantum fluctuations. 

At higher temperature the fluctuation contributions of quantum and thermal nature 
grow to a singularity in the neighbourhood of the critical temperature. The assump- 
tion that the Neel mean-field contributes for a major part to the magnetization and the 
susceptibility is no longer valid. Approaching the critical temperature T c the mean-field 
contribution to the magnetization goes to zero and strong diverging fluctuations are gen- 
erated at the one-loop order. This behaviour is common to the Heisenberg and XXZ 
magnetization. In addition the Neel order breaks SU(2) symmetry of the Heisenberg 
Hamiltonian inducing low momentum fluctuations near T c which is not the case in the 
XXZ-model. 

The influence of fluctuations decreases with the increase of the dimension D of the 
system due to the expected fact that the mean-field contribution increases relatively to 
the loop contribution. 

In dimension D = 2 the magnetization verifies the Mermin and Wagner theorem [EE] 
for 0, the fluctuations are larger than the mean-field contribution for any temperature. 
In a more realistic description another mean-field ansatz may be necessary in order to 
describe the correct physics. Indeed Ghaemi and Senthil [2.8. have shown with the help of 
a specific model that a second order phase transition from a Neel mean-field to an ASL 
(algebraic spin liquid) may be at work depending on the strength of interaction parameter 
which enter the Hamiltonian of the system. This confirms that another mean-field solution 
like ASL may be a better starting point than a Neel state when the temperature T 
increases. 
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Chapter 5. Mean-field ansatz for the 2d Heisenberg model 



We consider ordered quantum spin systems at finite temperature in which each lattice 
site is occupied by one electron with a given spin. Such a configuration can be constructed 
by means of constraints imposed through the specific projection operation [HE] which fixes 
the occupation in a strict sense. The constraint can also be implemented on the average 
by means of a Lagrange multiplier procedure jSJ. It is the aim of the present chapter to 
confront these two approaches in the framework of Heisenberg-type models. 

The description of strongly interacting quantum spin systems at finite temperature 
generally goes through a saddle point procedure which is a zeroth order approximation 
of the partition function. The so generated mean-field solution is aimed to provide a 
qualitatively realistic approximation of the exact solution. 

However mean-field solutions are not unique. The implementation of a mean-field 
structure is for a large part subject to an educated guess which should rest on essen- 
tial properties of the considered system, in particular its symmetries. This generates 
a major difficulty A considerable amount of work on this point has been made and a 
huge litterature on the subject is available. In particular, systems which are described 
by Heisenberg-type models without frustration are seemingly well described by ferro- 
magnetic or antiferromagnetic ( AF) Neel states at temperature T = fTS] ITS] . It may 
however no longer be the case for many systems which are of low-dimensionality (d < 2) 
and (or) frustrated [H21 EE1 EI] ■ These systems show specific features. An extensive anal- 
ysis and discussion in space dimension d = 2 has recently been presented by Wen [8U] . 
The competition between AF and chiral spin state order has been the object of very re- 
cent investigations in the framework of continuum quantum field approaches at T = 
temperature, see [78*1 

The reason for the specific behaviour of low- dimensional systems may qualitatively be 
related to the fact that low- dimensionality induces strong quantum fluctuations, hence 
disorder which destroys the AF order. This motivates a transcription of the Hamiltonian 
in terms of composite operators which we call "diffusons" and "cooperons" below, with 
the hope that the actual symmetries different from those which are generated by AF order 
are better taken into account at the mean-field level [SJ. 

In the present work we aim to work out a strict versus average treatment of the 
site-occupation constraint on systems governed by Heisenberg-type Hamiltonians at the 
mean- field level, for different types of order. This original work was published in [22] • 

The outline of the chapter is the following. Section l5~TI is devoted to the confrontation 
of the magnetization obtained through this procedure with the result obtained by means of 
an average projection procedure in the framework of the mean-field approach characterized 
by a Neel state. The same confrontation is done in section 15.21 for the order parameter 
which characterizes the system when its Hamiltonian is written in terms of so called 
Abrikosov fermions [SHI IE] i n d = 2 space dimensions. In section 15.31 we show that the 
rigorous projection [66J is no longer applicable when the Hamiltonian is written in terms 
of composite "cooperon" operators. 
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5.1 Antiferromagnetic mean-field ansatz 
5.1.1 Exact occupation procedure 

In this subsection we repeat the main steps of the partition function derivation given 
in section (|4.3|) with a slight modification. Instead of integrating over the Grassmann 
variables we construct a mean-field Hamiltonian "Hmf expressed in terms of the creation 
and annihilation fermion operators (p and /) and giving the same mean-field partition 
function Zmf as if one integrates over the Grassmann variables. Another point is the fact 
that instead of using the modified Matsubara frequencies as shown in subsection 13.3.41 we 
keep the imaginary chemical potential as an explicit term outside of the time derivation 

Starting with the nearest-neighbour Heisenberg Hamiltonian defined in equation (|4.3j) 



<i,j> i 

the partition function Z can be written in the form ()3.39j) 



Z {PFP) = I 

-4, CT (/3)=-&, CT (0) 

X>^ e -sm,<r&A) (5.2) 

where the {£,* a , £,i, a } ar e Grassmann variables corresponding to the operators {fj a , fi a } 
defined in section |3~T1 These Grassmann variables depend on the imaginary time r in the 
interval [0,(3]. The action S is given by 



SiitZrtZu,}) = [ P drJ2^Ar)d T ^(r) + W (PFP) ({^(r),^(r)})) (5.3) 
Jo 



where 



h {pfp) (r) =h(t)- hN(t) (5.4) 

and N(t) = EC(7( T )6,ff( r ) is the particle number operator. \x is the imaginary chemical 

potential introduced in section 13.21 describing the Popov and Fedotov procedure. The 
Hubbard-Stratonovich (HS) transformation defined in subsection 14.3.11 which generates 
the vector fields {&} leads to the partition function Z which can be written in the form 



Z (PFP) = J V ^^ . } y~S§dr { ^^^ (55) 
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In equation ()5.5|) the expression of Z is quadratic in the Grassmann variables {4* ff ,(;, ff } 
over which the expression can be integrated. The remaining expression depends on the 
fields {(pi{r)}. A saddle point procedure decomposes y?j(r) into a mean- field contribution 
and a fluctuating term 

<Pi(r) = (0i - SjJ + 5(pi{r) (5.6) 
where 0i are the constant solutions of the self-consistent equation 



as was shown in section l4~H for equation (|4.4fij) . 
The partition function takes the form 



Z(PFP) = z{ PF P) J V({m)e - 



S eff ({S&}) 



(5. 



where the first term on the right hand side corresponds to the mean-field contribution and 
the second term describes the contributions of the quantum and thermal fluctuations. 

In the following we focus our attention only to the mean-field part of the partition 
function 



Z (PFP) _ e - Seff , 
MF — e 



Tr 



> DrL MF 



(5.9) 



where Ti.^p P ^ is the mean-field part of the Hamiltonian ()5.4|) and reads 



where spin operators Si are given by (|3.2jl . (J -1 )^- is the inverse of the coupling matrix 
defined by (|4.2j) and N = Ylf'Lfiv * s the number operator. The mean-field Hamiltonian 

f|5.10|) expressed in terms of the creation and annihilation operators {f} a , fi, a } reads 



fi 



(PFP) 
MF 



= -^E( J_ %(^-^)-(^"^) 
+ E(4 /i 



2 



<£i_ 
2 



(5.11) 
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The mean-field Hamiltonian (|5.1H1 can be diagonalized by means of a Bogoliubov 
transformation as shown in appendix [01 and leads to 



+ £ + <5%-Ac-> } (5-12) 

i 

Fermion creation and annihilation operators {^r±\, A,(±)} are linear combinations of op- 



erators {/j T i(T , and the corresponding excitation energies uf^p reads 



{PFP) , K^i /r -I o\ 

The partition is then easily worked out and reads 

= e i ^ I(ft - 4)( *- S ' ) noosh(«) (5.14) 

i 

and the free energy is given by the expression 



(5.15) 



Going through the same steps as in section 14.5.11 the local mean-field magnetization 
[fhi j is obtained from 



rrii = 



and is related to the {<^}'s by 



rrii — tanh 



(PFP) 

(5.16) 



2(fi 



(5.17) 



and by virtue of the relations ()4.55|) and (|4.5fi|) one gets the self-consistent equation for 
the {rfii} 
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rhi = - y^Jfi 1 tanlT 1 (2m,) (5.18) 

3 n Jll 

If the local fields {-Bj} are oriented along a fixed direction e 2 , = mje^, the magnetiza- 
tions are the solutions of the self-consistent equations 

m i = -tanh ] (5.19) 



5.1.2 Lagrange multiplier approximation 

In section EIH we introduced the projector ()3.18|) 



R t = I V\/^ fUa ~ 1] (5.20) 



which allows to fix to one the number of spin- 1/2 per lattice site. Similarly to the preceding 
case in which the Popov and Fedotov imaginary chemical potential was used one can 
introduce the one-particle site occupation by means of a Lagrange procedure. In order to 
do that one has to replace the Popov and Fedotov projector P = -^e l ^ N by ()3.18|) . The 
Hamiltonian 7i then reads 

n (x) = ^./, ; .S;..s ; + Y^BJi + x ^ ~ l ) ( 5 - 21 ) 

i,j i i 

where Aj is a variational parameter and {rii = J2fLfi°} are particle number operators. 

(7 

Following the same lines as in section In. 1. II with the help of a Hubbard- Stratonovich 
transformation and a Bogoliubov transformation as shown in appendix [01 the mean- field 
partition function Z^ F can be worked out 

Zg), = e -*gtf (*-4)(*-4>^jj (1 + e - K(+))(1 + e - KH) {522) 

i 

with 



<+) = A + ( 5 - 23 ) 



<-) = A - M (5.24) 
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The parameter A is fixed through a minimization of the free energy with respect to A 



(A) 



MF 



{A,=A} 



(5.25) 



The minimization shows that the extremum solution is obtained for A = and 



(A) 



MF 



(ft - a,) . {0, - s) - m 2 co S h (Bm 



(5.26) 



which is different from the expression of equation (|5.15|) by a factor 1/2 in the argument 
of the cosh term. 

The magnetization can be obtained in the same way as done in subsection 15.1.11 One 
obtains 



(A 0) _ _ J-gjjJj 

1 2 



\ 



(A=0) \ 
/ 



(5.27) 



which is again different from the expression obtained in the case of a rigorous projection, 
see equation (J5.19J) . 

The uniform solutions m\ PFP ^ = (-l) ! m' PFP ' and m\ X ~ ^ = (— l) l m( A= °) for Jij = 
J $ {^i — ?j ± ff) have been calculated by solving the selfconsistent equations 

ffejai ,— ,ap} 

(15.19)1 and 1)5.27)1 . The results are shown in figure f)5.1j) . It is seen that the treatment of 
the site-occupation affects sizably the quantitive behaviour of observables. In particular 
it shifts the location of the critical temperature T c by a factor 2. Such a strong effect has 
already been observed on the behaviour of the specific heat, see refs. 112] . 

5.2 Spin state mean-field ansatz in 2d 

In 2d space the Heisenberg Hamiltonian given by 



h = — j Si.Sj 



(5.28) 



<i,j> 



can be written in terms of composite non-local operators {£V/} (" diffusons" ) [S] defined 

as 



v v - ftifj.l + fljhl 
If the coupling strengths are fixed as 



(5.29) 
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Figure 5.1: Magnetization vs. reduced temperature t = T/\J\. Full line: exact site- 
occupation. Dashed line: average site-occupation. 



where ff is a lattice vector {01,02} in the Ox and Oy directions the Hamiltonian takes 
the form 



where i and j are nearest neighbour sites. 

The number operator products {r^n,,} in Eq. (|5.31|) are quartic in terms of creation and 
annihilation operators in Fock space. In principle the formal treatment of these terms 
requires the introduction of a mean-field procedure. One can however show that the 
presence of this term has no influence on the results obtained from the partition function. 
Indeed these terms lead to a constant quantity under the exact site-occupation constraint 
and hence are of no importance for the physics described by the Hamiltonian (|5.31|) . As a 
consequence we leave it out from the beginning as well as the contribution corresponding 
to the {rii} terms. 

5.2.1 Exact occupation procedure 

Starting with the Hamiltonian 




(5.30) 




(5.31) 



<ij> 




(5.32) 



<ij> 
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Figure 5.2: Plaquette (□) on a two dimensionnal spin lattice 



the partition function Z can be written in the form ()5.2j) and the Hamiltonian in the form 
(|5.4|) . A Hubbard-Stratonovich transformation on the corresponding functional integral 
partition function in which the action contains the occupation number operator as seen 
in equation (|5.Hjl eliminates the quartic contributions generated by equation (|5.29|) and 
introduces the mean- fields {A^}. The Hamiltonian takes then the form 



n(PFP) = 2 ^ + ^ + A ^ p j.j _ ^ (5 33) 

' ' <«j> <«i> 

The fields {Ajj} and their complex conjugates { A^} can be decomposed into a mean-field 
contribution and a fluctuation term 



Aij = A* + SAij (5.34) 

The field Afj F can be chosen as a complex quantity A*f F = \Afj F \e 1 ^ . 

The phase </>^^ is fixed in the following way. Consider a square plaquette □ = 

(i, i + e x ,i + e x + e y ,i + e y j where e x and e y are the unit vectors along the directions 

Ox and Oy starting from site i on the lattice as shown in figure (|5.2j) . On this plaquette 
we define 



= E ^ F ( 5 - 35 ) 

(ij)eD 

which is taken to be constant. If the gauge phase 4>f- F fluctuates in such a way that 
stays constant the average of will be equal to zero in agreement with Elitzur's 

theorem [26]. In order to guarantee the 577(2) invariance of the mean- field Hamiltonian 
along the plaquette we follow EH1 |3U IHOj and introduce 
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g l - 1 ( !)' ? jf f. = j?. _|_ 



e *■ 4 ( ^ , if fj 



r*i + eV 



(5.36) 



where e x and join the site % to its nearest neighbours j. Then the total flux through 
the fundamental plaquette is such that = 7r which guarantees that the SU (2) symmetry 
of the plaquette is respected 

At the mean-field level the partition function reads 



where 



Z 



(PFP) 
MF 



-PH 



(PFP) 
MF 



(5.37) 



n 



(PFP) 
MF 



(5.38) 



After a Fourier transformation the Hamiltonian ()5.38j) takes the form 



n 



(PFP) 
MF 



Af, 



v\ 



+ EE 44 J ^ (PFP) 



keSBZ 



f k+f, 



with 



j^(PFP) 



+ AcosjZ^ kx , ky -iAsm^z-f kx , ky+7T 
+iAsm^z-f kx!ky+7r -fi- Acos^z-f kxjky 



(5.39) 



(5.40) 



where A = |A MF |. The Spin Brillouin Zone (SBZ) covers half of the Brillouin Zone (see 
figure EOl in subsection 14.2. 2|l . The 7^'s are defined by 



7* 



(cos k x ai + cos k v a 2 ) 



(5.41) 



where z = 4 is the coordination and M the number of sites. a\ and a 2 are the lattice 
parameters in direction Ox and Oy. The lattice parameters are not important for our 
study. We renormalize the wave vector k by the relations k x = k x .a\ and k y = k y .ci2 as 
shown in appendix El The momenta {k} act in the first half Brillouin zone (spin Brillouin 
zone) . 

Performing a Bogolioubov transformation which diagonalizes the remaining expression 
(|5.40|) in Fourier space leads to 



U 



(PFP) 
MF 



MzA 2 
\J\ 



k.a 



,( PFP ) at 



UJ (+),k,*P J (+)A*P(+),k,v + ^(-),fc,! /9 (-),fc, (J ^(-),fc^ 



(5.42) 
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The transformation is worked out in appendix ID. 21 The eigenenergies and 
are given by 



u t)Z = _/i + 2A t cos2 (^) + cos2 ( k v)} l/2 (5-43) 



and similarly 



= -fi- 2A[cos 2 {k x ) + cos 2 {k y )} 1/2 (5.44) 

The partition function Zmf has the same structure as the corresponding partition 
function in equation (|5.14jl and the free energy is given by 

fSP = ^jf - (2oosh/?A £ j) (5.45) 



k.o 



with 



6^ = 2 [cos 2 (k x ) + cos 2 (k y )] 1/2 (5.46) 

Finally the variation of Tmf with respect to A leads to the self-consistent mean-field 
equation 



k,a \ / 

with A( PFP ) = zA/\J\ which fixes A. 

5.2.2 Lagrange multiplier approximation 

Similarly to equation (j5.21|) one may introduce a Lagrange constraint and write 



(5.47) 



2 

<ij> <ij> 

In the mean-field approximation one have 



A; = A (5.49) 

and after a Bogoliubov transformation as for equation ()5.42|) the mean-field Hamiltonian 

reads 
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(A) _ AfzA 2 ^ / (A) t (A) t 

MF ~ | j | + ^(+),fc, CT P (+)fe l( T P (+)fe^ + ^(-),fe, CT P (-)fc, CT P (-)fc^ 



(5.50) 



with the eigenenergies 



uW Uff = A + 2A[cos 2 (/y + cos 2 (fcj,)] 1/2 



(5.51) 



and similarly 



= A - 2A[cos 2 (^) + cos 2 (^)] 



1/2 



(5.52) 



The expression of the free energy is now given by 



k,(i 



o) 



-0uj 

1 + e (+) 



(A) 



-flu 

1 + e (-) 



(5.53) 



As was shown in subsection 15.1.21 the minimization of this expression in terms of A 
delivers the solution A = and 



r 



(A) 



MF 



AfzA 2 
\J\ 



k,cr 



(5.54) 



where eg is given by equation (|5.46jl . 

The variation of T^j F with respect to A leads to the self-consistent mean-field equation 



>|J|e s AW' 



22 



(5.55) 



k.cr 



with AW = *A/|J|. 

Expressions in equation (|5.54j) and equation ()5.55|) should be compared to the expres- 
sions obtained in equation ()5.45|) and equation (|5.47|) . Figure IB~3l shows the behaviour 
of A for the two different treatments of site-occupation on the lattice. The exact oc- 
cupation procedure compared to Lagrange multiplier approximation doubles the critical 
temperature of the order parameter A as was shown in section 15.11 and in [Oj for the Neel 
state. 
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Figure 5.3: A vs. reduced temperature t = zT/\J\. Full line: exact site-occupation 
A( PFP \ Dashed line: average site ocupation A^ x \ 



5.3 Cooperon mean-field ansatz 

Starting from the Hamiltonian 

H = — \J\ ^ S{Sj 



(5.56) 



<i,3> 



one can introduce a further set of non-local composite operators {C^-} (" cooperons" ) 



fi,lfi,l - fi,ifj,1 



This leads to the expression 



riiTij 



(5.57) 



(5.58) 



where n« = Ylfiafi,v was explained in section IB~^1 for the Hamiltonian composed of 

a 

the operator T>^ the second term in the right hand side of equation ()5.58|) is quartic in 
the fermion creation and annihilation operator / hence one needs to use a Hubbard- 
Stratonovich transformation in order to reduce this term to a quadratic form. Calculations 
show that the second term is irrelevant. The solution of the self-consistent equations of 
the auxiliary field introduced by the Hubbard- Stratonovich procedure lead to {nj = 1}. 

5.3.1 Exact occupation procedure 

As in the preceding cases it is possible to implement a Hubbard-Stratonovich procedure 
on the {Cij} in such a way that the expresssion of the corresponding partition function 
gets quadratic in the fields {/»,-[, The corresponding HS fields are {Ty} and 
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1 1 <hJ> <h3> 



Introducing the homogeneous mean-fields {r^ = T = rj, one gets in Fourier space 



nj{PFP) 
rL MF 



E 

keBZ,, 



zT Tk 



-k,—a 



J k,a J -k,-a 



(5.60) 



with 



7fc = ^(cosk x + cosky) (5.61) 

The second term in this expression is complex since fi = in/ 2(3. In this representation 
it is not possible to find a unitary Bogolioubov transformation which diagonalizes H as 
shown in appendix ID .31 Hence a rigorous implementation of the constraint on the particle 
number per site is not possible. 

Reasons for this situation are the fact that the Hamiltonian contains terms with two 
particles with opposite spin created or annihilated on the same site which is incompatible 
with the fact that such configurations are not allowed in the present scheme. Terms of 
this type are typical in mean-field pairing Hamiltonians which lead to a non-conservation 
of the number of particles of the system as for the .BCS-normal transition. 

5.3.2 Lagrange multiplier approximation 

If the sites are occupied by one electron in the average the Lagrange procedure works 
opposite to the exact procedure. Here a Bogoliubov transformation can be defined and 
used to diagonalize the mean-field Hamiltonian. We do not develop the derivation of the 
mean-field physical behaviour here since it has been done elsewhere [3 El IS] • 



5.4 Summary and conclusions 
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5.4 Summary and conclusions 

In summary we have shown that a strict constraint on the site-occupation of a lattice 
quantum spin system described by Heisenberg-type models shows a sizable quantitative 
different localization of the critical temperature when compared with the outcome of an 
average occupation constraint. Consequently it generates sizable effects on the behaviour 
of order parameters, see also reference [22]. With exact site-occupation the transition 
temperature of antiferromagnetic Neel and spin states order parameters are twice as large 
as the critical temperature one gets from an average Lagrange multiplier method. 

Opposite to the average procedure the exact occupation procedure can not be used on 
cooperon states mean- field Hamiltonian. Cooperon are BCS pairs then destroy two quasi- 
particles in favor to create a new one which is a pair and then the number of particle are 
not concerved opposite to the exact occupation method. No fluctuations of the number 
of particles are tolerated by the exact occupation procedure. 

Due to the complexity of quantum spin systems the choice of a physically meaningful 
mean-field may depend on the coupling strengths of the model which describes the sys- 
tems |81j . A specific mean-field solution may even be a naive way to fix the "classical" 
contribution to the partition function which may in fact contain a mixture of different 
types of states. As already mentioned many efforts have been and are done in order to 
analyze and overcome these problems by means of different arguments I78[ l80j. 

In a more realistic analysis one should of course take care of the contributions of quan- 
tum fluctuations which may be of overwhelming importance particularly in the vicinity 
of critical points. Chapter El study the implications of the phase flutuations around the 
7r-flux mean-field state. 
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A Neel ansatz is not necessarily a good candidate for the description of two dimensional 
quantum spin systems. We showed in chapter 0] that such an ansatz breaks the SU(2) spin 
symmetry generating Goldstone bosons which destroy the Neel order. A better candidate 
seems to be the spin liquid ansatz since it keeps SU{2) symmetry unaffected. For this 
reason in the following we concentrate us on the spin liquid phase. 

Quantum Electrodynamics QED^+i) is a common framework which can be used to 
describe strongly correlated systems such as quantum spin systems in 1 time and 2 space 
dimensions, as well as related specific phenomena like high-T c superconductivity fI7\ EHl 
EOJ EZ|- A gauge field formulation of antiferromagnetic Heisenberg models in d — 2 
space dimensions maps the initial action onto a QED3 action for spinons |2H1 EZI- This 
description raises the problem of the mean-field solution and the correlated question of 
the confinement of test charges which may lead to the impossibility to determine the 
quantum fluctuation contributions through a loop expansion in this approach jH2l IHH1 E2] • 

We consider here the 7r-flux state approach introduced by Affleck and Marston El] • 
The occupation of sites of the system by a single particle is generally introduced by means 
of a Lagrange multiplier procedure [HJ Ej- In the present work we implement the strict 
site-occupation by means of constraints imposed through a specific projection operator 
which introduces the imaginary chemical potential proposed by Popov and Fedotov [HE| 
for SU(2) and modifies the Matsubara frequencies as explained in chapter El 

Here we concentrate on the generation and behaviour of spinon mass which stems 
from the presence of a U{1) gauge field. Appelquist et al. [SJ E] have shown that at 
zero temperature the originally massless fermion can acquire a dynamically generated 
mass when the number N of fermion flavors is lower than the critical value N c = 32 /V 2 . 
Later Maris [52J confirmed the existence of a critical value N c ~ 3.3 below which the 
dynamical mass can be generated. Since we consider only spin-1/2 systems, N = 2 and 
hence N < N c . 

At finite temperature Dorey and Mavromatos [21] and Lee 0H] showed that the dy- 
namically generated mass vanishes at a temperature T larger than the critical one T c . 

We shall show below that the imaginary chemical potential introduced by Popov and 
Fedotov [66J modifies noticeably the effective potential between two charged particles and 
doubles the dynamical mass transition temperature, in agreement with former work at 
the same mean- field level [22] • This original work published in [23] leads to another work 
which is not directly related to the PFP and concerns the introduction of a Chern-Simons 
term in our theory [20] . 

The outline of the chapter is the following. In section 16.11 we give a definition of a 
spinon. In section E21 we derive the Lagrangian which couples a spinon field to a U(l) 
gauge field. Section |p r 73l is devoted to the comparison of the effective potential constructed 
with and without strict occupation constraint. In section KTM we present the calculation 
of the mass term using the Schwinger- Dyson equation of the spinon. Section IfTH gives an 
outlook on the spin-charge separation phenomenon and confinement problems. 
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6.1 Definition of a Spinon 

Electrons are particles with charge e and spin 5 = 1/2. Formally they can be defined by 
creation and annihilation operators Cl a and CV j(J at position r and spin projection a = 
±1/2. In this description the t — J model which was largely studied for low-dimensional 
systems after the discovery of High-T c superconductivity reads [U EU EH1 E2] 

H = -YSj« C l C ^ + Y^JiA-Sj (6.1) 

ij a ij 

where tij is the so called hopping energy when the electrons jumps from j to site i and 
Jij is the spin coupling matrix as defined in section Kill 

It is believed that for a linear chain the t — J model leads to the phenomenon of spin- 
charge separation for strong electronic correlations [HU EH] • Observation of spin-charge 
separation in one-dimensional materials such as SrCu02 confirms this belief [HI]. In this 
framework the creation and annihilation operators C\ a and CV jCr are represented in terms 

of spinon ft a , f^ a , holon hi, h? and doublon dl , creation and annihilation operators 
and read jHH] 

CW = f^hl+e^f^d? (6.2) 

where f fulfils fermion anticommutation, h and d fulfil commutation relations. e> is the 
antisymmetric tensor. The transformation (jfi.2j) is exact when the occupation lattice sites 
of these entities is controlled by a constraint and reads 

hplr + fl^frA + fljrM + d^r = 1 (6.3) 

corresponding to the fact that each lattice site can be occupied by a hole, a spin a =|, a 
spin a =1 or double occupancy. This description is the so called slave-boson method. One 
sees that holons and doublons are charge excitations while spinons are spin excitations of 
the electonic systems. 

If the Popov and Fedotov procedure (PFP) is adopted holons and doublons are not 
allowed to live on the lattice and the constraint ()6.3|) reduces to 

4/fiT + fljrM = 1 (6.4) 

In this chapter since one considers the PFP only spinon entities are under consideration 
in our Heisenberg models describing cuprates in their insulating phase as explained in 
section 13.11 
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6.2 The 7r-flux Dirac action of spinons 

We have seen in section 15.21 that the Heisenberg Hamiltonian H = —J Si.Sj can 



<i,i> 



be expressed in terms of the diffuson operator {X*^} and leads to the 7r-flux mean-field 
Hamiltonian 



{PFP) 
MF 



n 

+ EE 

k&SBZ a 



Af; 



A 2 

'JT\ 



J k,a J k+if,<r 



-fi + A cos I z'j^^ -zA sin f^.^+TT 
+iA sin \ z-f kxtky+7T -fi-A cos \z~f kxiky 



■> k,a 



(6.5) 



for which the eigenvalues read ^^Pp = i 2Aa/cos 2 (A; x ) + cos 2 (k y ) and are shown in 
figure 16.11 





(a) (b) 
Figure 6.1: (a) Representation of the energy spectrum 



(PFP) 
(-)fc,<r 



— 2Ay / cos 2 (fc x ) + cos 2 (k y ) for k x and k y belonging to [— f,f] and (b) the contour 
representation of the energy spectrum showing the presence of the nodal points (±f , ±f ) 
where the energy is equal to zero. 

We are interested in the low energy behaviour of the quantum spin system described 
by the Hamiltonian (|6.5|) in the neighbourhood of the nodal points (k x = ±|, ±|) where 

the energy gap {^^P a — u P)pPj vanishes as shown in figure 16.11 

As already shown in earlier work by Ghaemi and Senthil [2H] and Morinari [2Zj the 
spin liquid Hamiltonian ()6.5j) for spin systems at low energy can be described by four- 
component Dirac spinons and the corresponding Dirac Hamiltonian reads 



keSBZ 



Ti 
T 2 



Ak_ 



T 2 
Ti 



^ka 



(6.6) 



k+ = k x + k y and k_ = k x — k y , A = 2A cos | and 
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( f la,kcr \ 
f 

flak 
V ^2bka / 



'lb.kcr 



(6.7) 



Indices 1 and 2 in the definition of ip given by equation ()6.7|) refers to the two independent 
nodal points f|, |) and (— f , f), indices a and 6 originate from a linear transformation 
between and /g +7f as detailed in appendix IE. 11 

The Dirac action of a spin liquid in (2+1) dimensions is derived in appendix IE. II In 
Euclidean space the action reads 



S E = dr d 2 f^ra W (d T -fi) + ^ k d k 



(6.8) 



where A = 2A cos | is the "light velocity" space component of the three "light velocity" 
vector v M = ( 1, A, A ) , and the {7 M }'s are the Dirac gamma matrices in (2+1) dimensions. 



6.2.1 "Gravitational" effects 

In the action (|6.8|) the "light velocity" which is affected by the parameter A can be 
seen as a space-time curvature parameter. This term must modify the covariant derivative 
due to the "gravitational" effect induced by the "light velocity" term. We show now in 
which way one must modify the theory in order to include this gravitational effect which 
finally is no longer important in our spinon action. 

In a curved space-time the Dirac action of a relativistic fermion reads [HE] 



S Di rac = J d z xE^e» ( <9„ + -uj^ ab [ 7 a , 7 b ] J V> (6.9) 

where are triads (also called dreibein) and E = det relates the "flat" space (with 
"flat" index m, where there is no gravitational effect) and the curved space (with "curvy" 
index /z, where there is a gravitational effect). The second term in brackets (equation 
(16. 9|) , where w^^b is the connection coefficient) comes from the preservation of the local 
invariance under Lorentz transformations in the presence of curved space-time jHHl (it is 
similar to the gauge invariance transformation of the covariant derivative — > D^). 

Comparing equation ()6.8j) with ()6.9)1 one shifts the imaginary time derivation d T — > 
d T + fx in ()6.8)1 which leads to a new definition of the Matsubara frequencies only for the 
fermion fields ip [66 . These frequencies read 



2tt 



= uf,u -n/i= -p(n + I/ 4 ) 



(6.10) 
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Further we define the dreibein as 
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Va 
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J 



Hence spinons move in a "gravitational" field which can be characterized by the metric 



T 

A 
A 



(6.11) 



One should add the term 



(6.12) 



in equation (J6.8() in order to verify the local invariance of the Lorentz transformation 
under the "gravitational" field. The metric can be handled jTHj assuming A = 1 without 
altering the physics of the problem as we will explain below. 

Since the metric (16.11(1 is no longer that of a flat space-time the connection coefficient 
Lo^ ab is defined by 



u^ab = - (df,e ab - T^e^b) (6.13) 

The Ps are Christoffel symbols (r^ 7 = \ (g^ + g m%v - g UJ ^) with = d^g^) and 
are the dreibein. Since A can be considered as constant in space-time we see clearly 
that the dreibein are also constant with respect to the space-time coordinates. Hence 
^a,ab = in a dilated flat space-time with the Euclidean metric (|6.11() . Since we showed 
that ()6.12j) is equal to zero one can put A = 1 in equation ()6.8j) and finally equation ((6.8(1 
reduces to 



Se = f Q dT J d2 ?Yj^™ (t) 7Afe ( r ) ( 6 - 14 ) 

Recall that we shifted the time derivative d T by the imaginary chemical potential \i mod- 
ifying the Mastubara frequencies (joMOjh More details concerning the modification of 
the Matsubara frequencies are given in subsection 13.3.41 This modification will induce 
substantial consequences as it will be shown in the following. 
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6.2.2 Quantum Electrodynamic Spinon action in (2+1) dimen- 
sions 

Since the Heisenberg Hamiltonian (j6.6|) is gauge invariant in the U(l) transformation 
ip — > e l9e ip, where 9 is a space-dependent scalar function, the Dirac action can be written 
in the form 







1 

— ( 

2 



S E = I dr I d 2 f <{ --a„ (f, r) [(□<$"" + (1 - A)^)] a, (f, r) 



+ ( r ) ^ _ ^ ( r ) r ( ^ 6 ' 15 ' ) 



Here (7 comes as the coupling strength between the gauge field a M and the Dirac spinons 
ip. The gauge field a M is related to the phase 9 (r) of the spinon at site r through the 
gauge transformation ff >a — > e tge ^f? i(7 which keeps the Heisenberg Hamiltonian invariant. 
From the definition of ip one gets ^Vo- — » e ige ^^p^ a . It is clear that # (r) is the phase at 
the lattice site r and that a M (r) = 0^9 (f). Hence the fluctuations of the flux (ft through 
the plaquette (see figure E2J are directly related to the circulation of the gauge field 
around the plaquette 

= g v in) - e (&)) 

<i,j>eO 
dl.a 

□ 

Phase fluctuations of lattice sites are measured by the gauge field a M . In (I6.15jl the first 
term corresponds to the "Maxwell" term —\}\ a ,f^ v of the gauge field a M where f^ v = 
d^a u — d v a^ A is the parameter of the Faddeev- Popov gauge fixing term —A (d^a^) 2 [HQ] , 
5^ u the Kronecker 5 and □ = d^. + V 2 is the Laplacian in Euclidean space-time. 

Since the spinons ip are minimally coupled to the gauge field a M , materializing the 
effect of fluctuations around the mean-field spin liquid ansatz described by the "Diffuson" , 
the energy-momentum conservation leads to consider the gauge-invariant and symmetric 
energy-momemtum tensor 



dC 

®^ = ^7^ — \ d " a s ~~ 9i*v£ + d P ifiipO-v) (6-16) 

where C is the QED 3 Lagrangian of the spinon deriving from Se = / drC in the 
Euclidean space. The Maxwell Lagrangian —\f liv f liv = (D5^ + (1 — X)d^d u ) a u is 
added to ()6.14j) in order to verify the energy-momentum conservation [HH] 



(6.17) 
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Notice that "photon" described by the gauge field are not present in nature and thus 
the dynamical Maxwell term above looks inadequate as it stands. In order to be coherent 
one should do the change of variable 



a^, = a^/g (6.18) 

where g is the coupling constant between the spinon and the gauge field. The spinon 
action ()6.15|) reads 



S E = j^dr J d 2 f I -±- 2 a, (f, r) [(W + (1 - \)d»d v )) a v (f, r) 

+ ^ypf a (r) [7^ (fy - ictp)] ipra (r) J (6.19) 

The dynamical term can now be removed in the limit of g — > oo verifying the energy- 
momentum conservation law ()6.17|) . In this limit we are treating a strongly correlated 
electron system. Keeping this in mind we go further in the development using ()6.15|) . 

6.3 The "Photon" propagator at finite temperature 

Integrating over the fermion fields tp, using relation (jA.24j) in appendix El the partition 
function Z a] = J V (ip, a) e~ SE with action Se given by equation ()6.15|) leads to the 
pure gauge partition function 

Z[a] = J V{a)e- S *ff [a] (6.20) 
where the effective pure gauge field action S e ff [a] comes in the form 



S, 



eff 



[a} = J\t j rf 2 f|-ia M [(D^ + (1-A)^9")]a l/ J- In det ^(^-i^)] 



(6.21) 



Here one recognizes similarity between the second term in (|6.21j) and the log-det term in 
()4.39j) . Following the same steps as in section 14.41 one can develop the second term in the 
effective gauge field action S e ff [a] into a series and write 



■o ^ 

In det [ 7 „ (0„ - iga^)} = In det Gp 1 - J2~ Tr [^7%]" (6.22 

1 ^ 

n=l 
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where G F l (k — k') = z ^^ (5(/c — k') is the fermion Green function in the Fourier space- 
time with k = (u)F,ni kj , hence Gp = — i'^ !L (2k) 2 (35 (k — k') . The first term on the right 

hand side of equation ()6.22|) being independent of the gauge field {a^} can be removed 
from the series since we focus our attention on pure gauge field terms. The first term 
proportional to the gauge field n = 1 in the sum vanish since tvy^ = 0. Keeping only 
second order terms in the gauge field in order to treat gaussian fluctuations one gets the 
pure gauge action 



Sfh [«] 
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dr Id r< 



-a,i{W + {l-\)d»d v )}a u 
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(6.23) 



The second term in equation ()6.23|) is worked out in details in appendix IE. 21 and one gets 



o(2) 



2(3 
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d 2 q 
(2k) 
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A (o)- 



+ n^(g) a v (q) 



(6.24) 



A 



(o)- 

[IV 



- [(□5 M1/ + (1 — \)d^d u )] is the bare photon propagator. The detailed calculation 
of the polarization function U.^ is given in appendix lE.2l Equations ()E.23|) and (|E.24|) give 
the components of IL^. The finite-temperature dressed photon propagator in Euclidean 
space (imaginary time formulation) verifies the Dyson equation 



A-J = AjT' + IV (6-25) 
Finally the gauge effective action reads = — j^Yl I ^§2 a iJ.(~l)^u (l) a ^( ( l)- 

U! B 

6.3.1 Comparison of the Popov and Fedotov procedure with the 
Lagrange multiplier method 

Table l6~T1 compares the polarization function components IT, II2 and II3 obtained by 
means of the Lagrange multiplier approximation for which A = as shown in section 15.21 
with the PFP as computed in appendix IE. 21 Differences appear in the denominator term 
D(X, Y) and in the numerator of the integrant in II 2 where the cosine terms are replaced 
by sine terms. 

We now push the comparison further and show below how the PFP modifies the 
effective potential between two test particles and also affects the dynamically generated 
mass of the spinons. 
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Table 6.1: Comparison of polarization function components obtained from the Lagrange 
multiplier approximation and the Popov and Fedotov procedure 
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Lagrange multiplier method 
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J^dx \n2D(X,Y) 



D(X,Y) 



cosh (j3q\J x{\ — x)j + cos(27rxm) 



cosh (j3q^x{l — x)j + sin(27rxm) 



6.3.2 Covariant description of the polarization function 

One may believe that a system at finite temperature breaks Lorentz invariance since the 
frame described by the heat bath already selects out a specific Lorentz frame. However 
this is not true and one can formulate the statistical mechanics in a Lorentz covariant 
form [T§] . 

We consider the three dimensional Euclidean space. Define the proper 3-velocity 
of the heat bath. In the rest frame of the heat bath the three velocity has the form 

= (1, 0, 0) and the inverse temperature f3 characterizes the thermal property of the 
heat bath. 

Given the 3-velocity vector w M one can decompose any three vector into parallel and 
orthogonal components with respect to the proper velocity of the heat bath, the velocity 
u^. In particular the parallel and transversal components of the three momentum g M with 
respect to w M read 



(6.26) 



(6.27) 



Similarly one can decompose any vector and tensor into components which is parallel and 
transverse to a given momentum vector g M 
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Figure 6.2: The dressed photon propagator. Wavy lines correspond to the photon and 
solid loops to the fermion insertions 



(q.u) 

= u^-^q^ (6.28) 

V = <W-^ir ( 6 - 29 ) 

q A 

It is now easy to define second rank symmetric tensors T Mi , constructed at finite temper- 
ature from g M , and 5^ v which verifies q^T^ = 



V = ^-u"u v -^ (6.30) 

q 2 _ 

= ~pU^u u (6.31) 

C, v = 6^-^- (6.32) 

where and verifies the relation 



V + ^V = ( 6 - 33 ) 

Since one considers a spin system at finite temperature and "relativistic" covariance 
should be preserved the polarization function may be put in the general form ^H] 



IV = U A V + n B 5^ (6.34) 

and the Dyson equation ()6.25|) can now be expressed in a covariant form if one uses 
relation (jFT%E)l . 

6.3.3 Dressed "photon" propagator 

Inverting the Dyson equation ()6.25|) the dressed photon propagator is obtained by 
summation of the geometric series 

V = a{2 + A(°).(-n^).A^ + AS2.(_ n ^).Ag.(-n'o.A55 + ■ • • (6.35) 

Figure E21 shows the Feynman diagrammatic representation of the Dyson series (jo.35j) . 
The dressed photon propagator reads 



Chapter 6. Two-dimensional Heisenberg model and dynamical mass 
74 generation in QED3 at finite temperature 



A _ Afj,v B^ v q^q u 

(6.36) 



where II a and IL3 are related to IT^ by 



IL4 = nx + n 2 (6.37) 

n B = n 3 (6.38) 

The expressions of IT, II2 and 1I3 are explicitly worked out in appendix IE. 21 Here q° = 
^-m is the boson Matsubara frequency energy component of the photon three- vector q^ . 

The dressed photon propagator is clearly expressed in a covariant form since the basis 
tensor A^ u and are as shown in subsection 16.3.21 Remarkably the dressed photon 
propagator is composed of a longitudinal and a transverse part with respect to the photon 

momentum unlike the bare photon propagator A$ — 8^ — (1 — l/A)^ 1 when 



one takes the Landau gauge fixing condition A — > 00. Even is transverse to the photon 
momentum q^. In the Landau gauge one gets 



q 2 + n A q 2 + U B (q 2 ) 2 

showing clearly that the dressed photon propagator presents a longitudinal part with 
respect to the photon momentum since q^A^ v = q^B^ v = 0. 



6.3.4 Effective potential between test particles 

The questions to which one answers here concern the derivation of the effective potential 
interacting between two spinons and the impact of PFP site-occupation constraint on it. 

The spinon action (|6.15|) describes a gas of spinons and photons coupled together by a 
coupling constant g which can be interpreted as the charge of the spinon. Paying attention 
to the analogy with a plasma made of electrons and the spinon gas one can identify 
the interaction between two spinon like one would identify the interaction between two 
electrons. The term proportional to the product q a -qb where q a and q b are the charges of 
particles a and b correspond to the Coulomb interaction between two charge and reads 



V{R = \f a - f b \) = -27rg a g b ln (|f a - r 6 |) (6.40) 

in a two-dimensional system where particles are at position f a with charge q a and r b with 
charge q b . The corresponding action comes as 



^plasma 



^ l ~ <7a<ft 27rln (\f a - f b \) + const. ^q 2 a I 

L a=£b a ) 
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for a plasma two-dimensional system at temperature 1/(3. Define the spinon three-current 
jp = —ig^j4> (r r in^>a where ip is the Dirac spinor of the spinon defined in section IB~2l In a 

classical picture the time component of the current vector reads jo(x) = ^q a 5 (x — x a ). 

a 

Identifying the charge q with the component jo of the three- vector the Coulomb inter- 
action comes as 



q2 

V (R) = q a .q b p) . m ... (n ; In Z (6.42) 
dj (R, t) dj (0, r) 

where Z oc e~ Splasma is the partition function of the corresponding action S p i asma . Now 
using relation ()6.42|) with the partition function Z[i(>,a] — J T> (i/j, a) c~ Se with the spinon 
action Se given by equation (j6.15|) and rewritten as 



S E = / / 3 dr [ d 2 f \ -~a„ (f, r) [(OP* + (1 - \)d»d v )} a v (f, r) 



o 



+ X^ CT "fAipe* (r) + j M (f, r) a M (f, r) J 



(6.43) 



Considering the term proportional in j^, one finally gets from deriving the effective po- 
tential 

V (R) = q a q b / dr<a (R, r) a (0, r) > (6.44) 
Jo 

Here < ao (R, r) ao (0, r) > is to identify with the dressed photon propagator time com- 
ponents A 00 . The effective static potential V(R) between two test particles (spinons) of 
opposite chages q a = —q b = g at distance R is given by 



V(R) = -g 2 I drA 00 (r,R) 
Jo 



(6.45) 



After a Fourier transformation on A o the time components of the dressed propagator 
A M!/ in equation ()6.45|) one gets 

V(R) = -9 2 ^j (0 A oo(?° = 0, q)e^ (6.46) 



with 



Aoo (g° = 0, q) = — f ^ (6.47) 

g 2 + n 3 [q° = fm = 0j 
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T = 0.5 




10 20 30 40 50 



R 



Figure 6.3: Effective static potential with (fat line) and without (dotted line) the Popov- 
Fedotov imaginary chemical potential for the temperature T = {0.001, 0.05, 0.5}. 



The effective potential then reads 



V(R) = -f dqqJ (qR).—^J- (6.48) 

where Jo(qR) is the zero order Bessel function of the first kind. 



The polarization contribution n 3 (g° = 0, q) is equal to ^ dx In (2 cosh f3q^x(l — r 



when taking the Popov and Fedotov imaginary chemical potential into account. This has 
to be compared to the expression ^| j Q l dxln ^2 cosh ^q^x(l — x)^j when the Lagrange 
multiplier method for which A = is used [21], a detailed comparaison is given in table 

o 

For small momentum q — ► 0, Ih^m = 0) can be identified as a mass term (M {PFP) (/5)) 2 
and reads 



limn 3 (m = 0) = (M (PFP) (/3)) 2 = In 2 (6.49) 
q^0 Tip 



For R ^> (Mq PFP ^) 1 the effective potential reads 



vat,® K -t. r iq . ^ 



2tt 



f+(Mr p) y 



NV 8 nRMr P) 
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Figure 6.4: Temperature dependence of the dynamical mass generated with (full line) and 
without (dashed line) the use of the Popov-Fedotov procedure. 

where N = 2 since we consider only S — 1/2 spins. 

Figure l6~3l shows the effective potential between two opposite test charges at distance 
R > (M { Q PFP) )-\ The screening effect is smaller when the imaginary PFP chemical 
potential \x is implemented rather than the Lagrange multiplier A. By inspection one sees 
that (Mf^)" 1 = V2(M^y\ The main effect of combinating of the PFP and thermal 
fluctuations is to increase the range of effective static interaction, which is infinite at zero 
temperature, between two test particles of opposite charges compared to the case where 
the Lagrange multiplier method is used. 

6.4 Dynamical mass generation 

Appelquist et al. Ej have shown that at zero temperature the originally massless 
fermion can acquire a dynamically generated mass when the number N of fermion flavors 
is lower than the critical value N c = 32/7r 2 . Later Maris confirmed the existence 
of a critical value N c ~ 3.3 below which a dynamical mass can be generated. Since we 
consider only spin- 1/2 systems, N = 2 and hence N < iV c . 

At finite temperature Dorey and Mavromatos j21j and Lee jlH] have shown that the 
dynamically generated mass vanishes at a temperature T larger than the critical one T c . 

We now show how the Popov and Fedotov procedure doubles the "chiral" restoring 
transition temperature T c of the dynamical mass generation. The Schwinger- Dyson equa- 
tion for the spinon propagator at finite temperature reads 

G~\k) = G<»-\k) - J -^- 2l Mp)^(k - P)T„ (6.50) 

b>F,n 
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where p = (p = LOp tTl ,P), G is the spinon propagator, the spinon-" photon" vertex 
which will be approximated here by its bare value g^i v and A^ u the dressed photon propa- 
gator ()6.36j) . The second term in (|fi.5()|) is the fermion self-energy S, (G- 1 = Gty — E). 
Performing the trace over the 7 matrices in equation (|fi.5()j) leads to a self-consistent 
equation for the self-energy 



<+>F,n 

E (k) corresponds to a mass term which can be estimated at low energy and momentum 
limit m{(3) = E(/c) ~ E(0). In this regime equation ()6.51|) simplifies to 



Admitting that the main contribution to ()6.52j) comes from the longitudinal part Aoo(0, — P) 
of the photon propagator ()6.52|) goes over to 



ily 



d 2 P 



(2tt) 2 \p2 + n 3 (m = 0) 'Zj\ n + P 2 + m(/?) 2 



(6.53) 



Performing the summation over the fermion Matsubara frequencies uJF, n the self-consistent 
equation takes the form 



a 



AttN 



dP\ 



Ptanh(3JP 2 + m(/3) 2 



P 2 + n 3 (m = 0) \ P 2 + m((3) 2 



(6.54) 



Relation ()6.54|) can be solved numerically with a cutoff A to control the ultraviolet 
integration limit 



AixN 



dP 



P tanh p J P 2 + m(f3) 2 



P 2 + n 3 (m = 0) JP 2 + m(f3) 2 



(6.55) 



It has been shown elsewhere |3] and [25 that in the limit of A going to 00 numerical 
results are stable and finite. Figure 16.41 compares the temperature dependence of the 
dynamical mass generated with and without the imaginary chemical potential introduced 
by the Popov and Fedotov procedure where a/A = 10 5 . By inspection of equation ()6.54)1 
and the corresponding result obtained by Dorey and Mavromatos [21] and Lee ^H] one 
sees that the imaginary chemical potential used which fixes rigorously one spin per lattice 
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site of the original Hamiltonian (|5.31|) doubles the transition temperature. This result is 
coherent with the results obtained elsewhere [22] and shown in section I5~2l where spinons 
are massless. 

Since the mass can be related to the spinon energy gap and m (T = 0) corresponds 
to the spinon-antispinon condensate < tp^oa^k^ a > an amoun t of energy at least equal 
to m(T) is necessary in order to break a pair of spinon-antispinon and liberate a spinon. 
Interpreting the spinon as the spin excitation part breaking a Cooper pair one can identify 
m(T) with a superconducting gap and evaluate the reduced energy gap parameter r = 
2 j™j) where m(0) is the mass at zero temperature and T c the transition temperature 
for which the mass becomes zero. Dorey and Mavromatos [23] obtained r ~ 10 Lee 
[IB] computed the mass by taking into account the frequency dependence of the photon 
propagator and obtained r ~ 6. We have shown above that the imaginary chemical 
potential doubles the transition temperature so that the parameter r is ~ 4.8 for a/A = oo 
to be compared with the result of Dorey and Mavromatos and r ~ 3 to be compared with 
Lee's result. Recall that the BCS parameter r is roughly equal to 3.5 and the YBaCuO 
parameter r ~ 8 as given by the experiment [72j . 



6.4.1 Antiferromagnetic Neel order parameter 

There is another possible physical interpretation of the dynamical generated mass different 
from the former one [31 1231 HE] • One may also consider it as the antiferromagnetic Neel 
order parameter as will be explained belove [2*8} """""} I5T1 162j . 

Consider the mean-field Hamiltonian for which one takes into account both the diffuson 
()5.38)1 and the Neel ansatz (|5.1()j) . This mean- field Hamiltonian reads 



'■j 



(6.56) 



<ij> 



<ij> 



Following sections 16.21 and 15.11 and after some transformation steps on ()6.56)1 one gets the 
Dirac action [2Zj 



S E = dr SrYj^ra 7° {d T - v) + A 7 fc ^ 



4>f 



(6.57) 



The "mass" term m^ e \ oc is similar to the dynamical generated mass m(/3) introduced 
above but not equal to m(/3). The Neel ansatz introduces the Pauli matrix a affected to 
the mass term rn-Neei-o 2 of the spinons. In the context of the insulator compound described 
by the Hamiltonian ()5.31j) the "crural" symmetry breaking term mn&i.? corresponds to 
the development of Neel order |4"H as can be understood by simple comparison of 
(j6.15|) with ()6.57|) . One sees that adding the Neel anstaz, pictured by the term (pi.Si in 
()6.56|) . changes ()6.15|) into ()6.57|) thus mNeei is clearly related to the Neel order parameter. 
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As a consequence of this derivation the dynamical mass generation introduces natu- 
rally the concept of second order transition from a Neel phase to a genuine paramagnetic 
spin liquid in two dimensional quantum antiferromagnets contrary to general wisdom, as 
explained in [2Z1I2E!- This general wisdom claimed that a phase transition can only take 
place on one hand between collinear magnets and Valence Bond Solid (VBS) paramag- 
nets and on the other hand between non-collinear magnets and spin liquids. 

Finally one can interpret the dynamically generated mass m(/3) obtained in section 
16 .41 as the emergence of something like a Neel phase from a spinon gas describing the spin 
liquid state, more precisely from a 7r-flux state. One may push further and say that if 
this interpretation is correct then one confirms the results obtained in section I5~T1 and in 
[HI 122] concerning the doubling of the Neel order transition temperature T c . Indeed the 
dynamical generated mass sees its transition temperature doubled by the PFP compared 
to the Lagrange multiplier method like for the Neel order parameter as shown in chapter 

E 

6.5 The PFP and the confinement problem : outlook 

The QED 3 theory described above deals with noncompact "Maxwell" theory, in other 
words the integration over the gauge field a M goes within the limits ]— 00, oof and the 
Maxwell action can be described on a lattice by 



noncompact 




(6.58) 



X,[iU 

In a rigorous derivation on the lattice the "Maxwell" action is in the compact form [HH] 



5 comparf = ^(l-cos/^) (6.59) 

X,jjLV 

with f x ^ v = a x>IM + a x+li>v - a x+u ^ - a Xi „ and -it < a x ^ < 7T. Using this compact version 
of the Maxwell theory the partition function reads 



/ l[da x ^e — (6.60) 

and takes the periodicity of the action into account. n liF is an integer coming from the 
periodicity of the cosine and introduces new entities called instantons. Polyakov showed 
[53] that the compactness causes important changes in physical properties. Indeed a pure 
compact Maxwell theory confines test charged particles in (2+1) dimensions due to the 
formation of an gas of instantons. The electrostatic potential between two test particles 
behaves like V(R) oc R at zero temperature and (2+1) dimensions while the electrostatic 
potential is deconfining in the non-compact form and behaves like V(R) oc lni? at zero 
temperature and in (2+1) dimensions. 
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However when matter fields are taken into account the system may show a decon- 
fined phase even in the compact formulation of the theory One may ask now 

whether spinons deconfine or not ? A large amount of work has been devoted these 
last decades to the search of an answer [TBJ E2 EH ED EH ESI E2] • The possible existence 
of a confinement-deconfinement transition could leads to an explanation for spin-charge 
separation in strongly correlated electronic systems [SU]- I n the present context one may 
ask how the Popov and Fedotov procedure affects a compact QED 3 theory of spinons. As 
seen in chapters El -El previous work has shown that the PFP probably modifies the physi- 
cal behaviour in a quantitative way but does not produce a deep qualitative modification. 
These questions are potentially open for further work. 

Previous works showed that deconfined spinons (minimally coupled with a compact 
U(l) gauge field) can appear in the region of "a second order quantum phase transition 
from a collinear Neel phase to a paramagnetic spin liquid in two dimensional square lattices 
with quantum antiferromagnets and short ranged interactions" [2B1 US] • In this case the 
Hamiltonian is the one we used in our previous description with "diffusons" to which we 
add the term - ^2 Si.Sj where g is a parameter which controls the importance of 

<i,j>£N.N. 

the antiferromagnetic ordering. When g — > oo we reach the algebraic spin liquid phase 
and when g — > one gets the Neel phase. There is a critical coupling parameter g c under 
which the spinons are gapped, condense and are confined (through a Higgs mechanism). 
This corresponds to the Neel phase. For g larger than g c spinons are again confined but 
this time due to proliferation of instantons [ZH]. However these results are obtained for 
zero temperature and experiments have not yet proven the existence of spinons. 
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6.6 Summary and conclusions 

We mapped a Heisenberg 2d Hamiltonian describing an antiferromagnetic quantum spin 
system onto a QED^+i) Lagrangian coupling a Dirac spinon field with a U(l) gauge 
field. In this framework we showed that the implementation of the constraint which fixes 
rigorously the site occupation in a quantum spin system described by a 2d Heisenberg 
model leads to a substantial quantitative modification of the transition temperature at 
which the dynamically generated mass vanishes in the QED^+i) description. It modifies 
consequently the effective static potential which acts between two test particles of opposite 
charges [23 . 

The imaginary chemical potential [66J reduces the screening of this static potential 
between test fermions when compared to the potential obtained from standard QED^+i) 
calculations by Dorey and Mavromatos [21] who implicitly used a Lagrange multiplier 
procedure in order to fix the number of particles per lattice site EI] since A = at the 
mean-field level. 

We showed that the transition temperature to "chiral" symmetry restoration corre- 
sponding to the vanishing of the spinon mass m(/3) is doubled by the introduction of the 
Popov-Fedotov imaginary chemical potential. The trend is consistent with earlier results 
concerning the value of T c [22] • It reduces sizably the parameter r = ^ffi determined by 
Dorey and Mavromatos [21j and Lee [38j . 

Marston [53] showed that in order to remove "forbidden" £7(1) gauge configuration of 
the antiferromagnet Heisenberg model a Chern-Simons term should be naturally included 
in the QED3 action and fix the total flux through a plaquette. When the magnetic 
flux through a plaquette is fixed the system becomes 27r-invariant in the gauge field 
and instantons appear in the system. This is the case when the present non-compact 
formulation of QED 3 is replaced by its correct compact version [Ul] I65j. 

The implementation of a Chern-Simons term [23] in a non-compact formulation of the 
spinon system constrained by a rigorous site occupation has been submitted to publication 
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There is a general consensus on the phase diagram of high-T c superconductivity. The sim- 
plest model which seems to take account of the strong correlation physics of the high-T c 
superconductors is the Hubbard model and its strong coupling limit the t — J model. The 
schematic phase diagram of high-T c superconductors shows an insulating antiferromag- 
netic phase in the underdoped regime and is well described by Heisenberg-like models. 
A superconducting phase is present for low temperature and for a finite range of doping 
with holes (or electrons) as shown in figure 12.11 Anderson presented the viewpoint of 
the existence of a spin liquid state for a possible key to understand the physics of highly 
correlated superconducting phase [I]. The concept of Resonating Valence Bond (RVB) 
states was further developed [JBl EH] and gave rise to the notion of spin-charge separation. 
The inclusion of fluctuations around the corresponding "diffuson" mean-field led to a U(l) 
gauge theory [213 I2H EH EI] which can be treated at the one-loop level. 

One of the approaches to high-T c supercontuctivity phenomena consists in a slight 
doping of the superconductor materials starting from the undoped insulating antiferro- 
magnet phase to the underdoped regime. The underdoped superconducting should retain 
part of the correlation features of the insulating phase. For this reason it is of high interest 
to study two dimensional antiferromagnetic strongly correlated spin systems. 



In chapter |U we presented and discussed applications of the mean-field and loop ex- 
pansion to the determination of physical properties of antiferromagnetic Heisenberg-type 
systems (supraconductors in the antiferromagnet insulating phase) in spatial dimension 

Dim 

We worked out the expression of physical observables (magnetization and suscepti- 
bility) starting from a specific mean-field ansatz and including contributions up to first 
order in a loop expansion in order to investigate the effects of fluctuation corrections to 
mean-field contributions at the gaussian approximation. The mean-field was chosen as a 
Neel state which is an a priori reasonable choice for spin systems described in terms of 
unfrustrated bipartite Heisenberg model. The results were compared to those obtained 
in the framework of spin wave theory. 
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The number of particles per spin lattice site was fixed by means of the rigorous con- 
straint imposed by the imaginary chemical potential introduced by Popov and Fedotov 

At low temperature the magnetization and the magnetic susceptibility are close to the 
spin wave value as expected, also in agreement with former work [9 . 

At higher temperature the fluctuation contributions of quantum and thermal nature 
grow to a singularity in the neighbourhood of the critical temperature. In addition the 
Neel order breaks SU (2) symmetry of the Heisenberg Hamiltonian inducing low momen- 
tum fluctuations near T c which is not the case in the XXZ-model. 

The influence of fluctuations decreases with the dimension D of the system due to the 
expected fact that the mean-field contribution increases relatively to the loop contribution. 

In dimension D = 2 the magnetization verifies the Mermin and Wagner theorem jHSj 
for T^O, the fluctuations are larger than the mean-field contribution for any temperature. 
In a more realistic description another mean-field ansatz may be necessary in order to 
describe the correct physics. Ghaemi and Senthil [2E] introduced a specific model in 
which a second order phase transition from a Neel mean-field to an ASL (algebraic spin 
liquid) may be at work depending on the strength of interaction parameter which enters 
the Hamiltonian of the system. These considerations enforce the belief that another 
mean-field solution like ASL may be a better starting point than a Neel state when the 
temperature T increases. 

In chapter we worked out a rigorous versus average treatment of the occupation 
constraint on spin systems governed by Heisenberg-type Hamiltonians at the mean-field 
level, for different types of spin mean-field ansatz. 

We showed that a strict constraint on the site occupation of a lattice quantum spin 
system described by Heisenberg-type models shows a sizable quantitative different lo- 
calization of the critical temperature when compared with the outcome of an average 
occupation constraint. With an exact site-occupation by spin the transition temperature 
of antiferromagnetic Neel and spin liquid states order parameters are twice as large as the 
critical temperature one gets from an average Lagrange multiplier procedure. 

The exact occupation procedure cannot be applied to a so called cooperon state mean- 
field Hamiltonian. Cooperons are .BCS-like pairs of particles. In this scheme the number 
of particles is not conserved, hence it is incompatible with a strict site occupation con- 
straint. 

In a further step we mapped a Heisenberg 2c? Hamiltonian describing an antiferromag- 
netic quantum spin system onto a QEDr 2 +i) Lagrangian which couples a Dirac spinon 
field to a U{\) gauge field. We considered the 7r-flux state approach introduced by Af- 
fleck and Marston [21 EI]. We implemented the strict site-occupation by means of the 
constraint imposed through Popov and Fedotov's imaginary chemical potential [6H] for 
SU(2) and used the modified Matsubara frequencies. In this framework we showed that 
the implementation of the PFP constraint which fixes rigorously the site occupation in 
a quantum spin system described by a 2d Heisenberg model leads to a substantial quan- 
titative modification of the transition temperature at which the dynamically generated 
mass vanishes in the QEDm+x) description. It modifies consequently the effective static 
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potential which acts between two test spinons of opposite coupling charge g. 

The imaginary chemical potential 66J reduces the screening of the static potential 
between test fermions when compared to the potential obtained from standard QEDm+i) 
calculations by Dorey and Mavromatos |24j . 

We showed that the transition temperature corresponding to the vanishing of the 
spinon mass m{(3) is doubled by the introduction of the PFP. The trend is consistent 
with earlier results concerning the value of T c [2 2 3 . It reduces sizably the ratio r = 
between the superconducting gap m(0) and the transition temperature T c determined by 
Dorey and Mavromatos [2I| and Lee [4*8*] . 

We conclude from this work that the Popov and Fedotov procedure which treats ex- 
actly the constraint of strict site-occupation does not modify qualitatively but quantita- 
tively the physical results obtained by means of the average Lagrange multiplier method. 



The QED3 theory described in chapter [""] deals with a non-compact version of the 
Abelian gauge field action. In order to remove "forbidden" U(l) gauge configuration of 
the antiferromagnet Heisenberg model a Chern-Simons term should be naturally included 
in the QED 3 action "S3 [2D]. 

The compact Maxwell theory confines test charged particles in (2+1) dimensions due 
to the formation of an instanton gas [I""""]. When matter fields are taken into account the 
system may show a confinement/deconfinement transition [J""2|. The possible existence of 
a confinement-deconfinement transition could give an explanation for spin-charge separa- 
tion in strongly correlated electronic systems [60J. The fundamental question concerning 
confinement invalidating a loop expansion is up to now unsettled ["""*"] l"""2~] . 

It would be interesting to implement a Chern-Simons term "2*3 EH] in a compact U(l) 
gauge formulation of the spinon system constrained by a strict site-occupation imposed 
by the Popov and Fedotov procedure. 
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Appendix A 

Grassmann algebra and coherent 
states 

Here we review the Grassmann algebra and some properties of coherent states. 

• The anticommutation relation of the fermion creation and annihilation operators 
are 

{fa,fl}=6 afi (A.l) 

{/«,/*} = (A.2) 
{flfl}=0 (A.3) 

• The anticommutation relations of the Grassmann variables with themselves and the 
creation and annihilation fermionic operators are 



{UZx}=Ux + U a = (A.4) 

= (A.5) 

= (A.6) 

{^CaT = A is a complex number (A. 7) 

{£«, fe} = {U, 4) = {CM = {C4) = o (A.8) 

where £* is the conjugate of £ a . Consequently the product of two identical Grass- 
mann variables is zero due to the anticommutation property 



e a = (Cj = (A.9) 

Conjugation rules : 
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&r = c (a. io) 

(CY = £a (A.ll) 



Hermit ic conjugation : 



(£/) f = fC (A.12) 



Derivation : 



!«*« = |(-«-: 



-r (a. 13) 



Define the operator O 



o(e,e*) = a+/3e+ 7 r+Are (a.w) 

where a, /?, 7, A are complex number. With this definition 



JW,0 = (A.15) 

= 7+Ae (A.16) 



d d 



° d :0{^C) (A.17) 



d?d£ 



Integration : 



d£l = j dC 1 = (A.18) 
dU = [dCC = l (A. 19) 



Using the definition (jA.14|) of the operator O 
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J dCO(C,C) = /3-At (A.20) 

J dCO(^C) = 7 + A (A.21) 

J d?dZO(t,C) = -A 

= - J «*0(£,f) (A.22) 

Notice that the integration operator is equivalent to an ordinary derivation operator. 

Gaussian integrals : 
For commuting variables 

/ IT o e °' A ° = detMpe^ (A.23) 

y - LJ - 2m 



For Grassmann variables 



/ JJ t^Xae Q ' A Q = [detM]e"- A ' (A.24) 

Coherent states : 

The fermionic coherent states |£ > are defined as 



|f> = e _E «^|0> (A.25) 
<f| = <o|e E « r/a (A.26) 



Application of the creation and annihilation operator on coherent states leads to 
the following expressions 
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/a|C> = Ut> (A.27) 
/llO = flH(l-Wl)\0> 

A 

= /i(i-a/i)n( i -^/i)io> 



A^a 



> 



9 



(i - ut) n (! - ^/i) io> 



at 

= (A ' 28) 

<ei/« = ^:<ei (A.29) 

<ei/J = <eiC (A.30) 

The overlap of two coherent states is 



= <0|n( 1 + G/«)(l-d/l)|0> (A.31) 

= e E «^' (A.32) 

and the elements of fermionic operators between two different coherent states |£ > 
and |£' > read 



< d\0 {{f a }, ie' >= e^'o ({£ Q }, {&}) 



(A.33) 



Appendix B 

Spin Brillouin Zone 

B.l Two dimensional bipartite lattices 

Here we construct the Spin Brillouin Zone (SPZ) of a two dimensional lattice. We de- 
fine the components of the wave vector k on the orthonormal reciprocal sublattice basis 
(ei, e.2 j which is given by eVe} = Sij where the direct basis is defined as 



ei 
e 2 



V2 



(B.l) 
(B.2) 



Since the direct basis is already orthonormal the reciprocal basis is identically equal to the 
first one ^e"i = e\,e2 = e^J. Figure IB~T1 shows a two dimensional bipartite lattice where 
blue points refers to one type A and red points to an other type B of sublattices. The 
crystal basis (e x , e y ) as well as the direct basis (ei, £2) are shown. 



■O 




•O 



O- 



•o 



Figure B.l: Two dimensional bipartite lattice with the lattice basis (e x , e y ) and the spin 
sublattice basis (ei,e*2) in direct space. 
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If one sets the lattice parameter a to 1 the wave vector k is defined by 

k = kiei + k 2 e 2 
h 



fci = 2tt 



2tt- 



A(B),1 
h 



(B.3) 
(B.4) 

(B.5) 



where N"a(b),1{2) is the number of sublattice type-A(.B) sites in the direction 1(2) with 

Na{B),1 A/"a(_b),1 



h e 



2 ' 



G 2Z and l 2 G 



AQl(B),2 Na(B),2 
2 ' 2 



G S forming the Spin Brillouin 



Zone depicted in figure POl The 7^ function defined in section 14.2.21 becomes 



Ik 



- e** 4 
v 

- cos k.i 



'ki + k 2 \ (k\-k 2 
cos 1 7= — I + cos 



V2 



V2 



(B.6) 



where k belongs to the Spin Brillouin Zone and D = 2. The total number of lattice sites 
is Af = Na,i-Na,2 + Nb,\-Nb,2- 

B.2 Three dimensional bipartite lattices 

In three dimensions the direct spin sublattice is introduced for the bipartite lattice. For 
the Neel state a face-centered cubic lattice is used. The lattice basis (e*i, e 2 , e^) reads 





1 






e\ = 


72 




+ e y 




1 






e 2 = 


7! 




+ e z 


e 3 = 


1 

71 


(e x 


+ e z 



(B.7) 
(B.8) 
(B.9) 



The reciprocal basis {e^} of the lattice is defined by e^e} = 5y and reads 
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These vectors are those of a body-centered cubic lattice. The wave vector k reads 



k = 
ki = 



hei + k 2 e 2 + k 3 e 3 
U 

2tt 



Na(B),i ^A{B),i 



where Ma(b\% is the number of site in the direction of the sublattice A(B). 
The 7g function for k in the three dimensional Spin Brillouin Zone reads 



(B.13) 
(B.14) 

(B.15) 



Ik 



- [cos (ki - k 2 + k 3 ) + cos (ki + k 2 - k 3 ) + cos (-h + k 2 + k 3 )] (B.16) 
o 



The total number of lattice is Af = YV^A,i + H^b,*- 




Figure B.2: Three dimensional bipartite lattice with direct space basis (ei,e 2 ,e 3 ) 

Figure IB. 21 shows a three dimensional bipartite lattice where blue points refer to one 
type of sublattice A and red points to the other type B. The crystal basis (e x , e y , e z ) as 
well as the direct basis (ej, e 2 , e 3 ) are indicated. 
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Appendix C 

Beyond the mean field : one-loop 
contributions 

In section 14.31 we denned the mean-field contribution of the matrix M as 



a 



Op,q 



det 



1 1 



dct 



Tq- [lUJ F>p + y Z i(uF,p ~ U F , q )] 6 p>q 



(C.l) 



where det G r 



is given by 



yj(^(p-g=Q)) 



and the contribution to the fluctuations 6<p 



M 



\5y~l(u) FyP - u F>q ) -\5(fl{u F>p - u F>q ) 



(C.2) 



It comes out that the term In det /3M can be decomposed into a series expansion 



ft °° 1 

lndet/3M = - V In 2 cosh ^||^(cj b = 0) II + Tr{V -{GqM^Y} 
z — ' 2 n 



(C.3) 



n=l 



where the fluctuation contributions appear in the second term of (|C.3|) . 



C.l First order contributions to the fluctuations and 
mean-field equation 

Here we aim to work out the n = 1 contribution to the sum of lndet/?M, Tr [GqMi], 
and extract the mean-field equation with respect to the Hubbard-Stratonovich auxiliary 
mean-field (p. The first order term in the fluctuations 5ip reads 
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[G M 1 



p,g 



Ai(p,q) Bi(p,q) 
Ci(p,q) Di(p,q) 



(C.4) 



where i stands for the position on the spin lattice and p and q refer to the fermion 
Matsubara frequencies up,p = ^-(p+ 1/4) and similarly for Zjp A . The matrix elements 
A,B,C and D are given by 



A(p,q) 



2 det G 



p i 



uf, p - y(^-(p-9 = 0)) 



<Wj(p-?)) 



+ 2deWf (j( P -^^ n 2 f(p-q)) 



Bi(p,q) 



Ci(p,q) 



Di(p,q) 



2 det G 



V L 



<ty>i ("j(p-?)) 



1 



2 det G, 



fyfi-jj-ip-q)) 



i 



2 det G„ 



¥(y(p-!)) 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



The first term in the sum of In det (3M reads 



Tr [GqMi] = 22[^(p,p) + A(p,p)] 



« p 



i p 



2detG 



P L 



#(0)*tf(0) + j (^(0)^-(0) + ft-(0)W(0)) 



2 det G r 



^(0)5^(0) 



(C.9) 



The value of ^X ^dctc J can be found with the help of the sum relation in [3U] and leads 



to 
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Tr [G Mi] = -§5>nh (f ||&(0)|| ) 



(CIO) 



The first term in order of the fluctuations Sip of the auxiliary field action So [<p\ in equation 
reads 



5 S [ip] 



Sip 



S<p=0 



_S_ 

Sip 



0* 



\\(pj - B^j S(u B = 0) + Sipj(u B ) 
f E J ^ [(& - ^) + " ^) ^(0)1 (C.ll) 



8ip=0 



Adding the first order term in the fluctuations Sip of Sq [ip] in (jC.10|) leads to 



SS, 



eff 



Sip 



8lp = 6 -^Ms^ 







Sip 

£ 



2 E^' (A - 4 



tanh | — 



.SiptiujB = 0) = 
(C.12) 



From ()C.12j) one obtains directly the mean-field equation of the Hubbard- Stratonovich 
equation ()4.4fij) . 



C.2 Second order fluctuation contributions 

The term n = 2 in equation (j4.43|) is given by 
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Tr[G M 1 G M 1 ] = ^[^^(wl + B^^p) 



l p,q 

+ Ci(p,q)Bi(q,p) + Di(p,q)Di(q,p) 
= [(^ + / 2 ) 2/C ^s S(Pi(-u B )6(Pi(u B 

i ui B 

+ \ [(<p~ /2) 2 K^ B ] 6<pt(-u> B )8<pt{wB) 



t,UJ B 



+ 1 [<Pi (fi + «<>b) ^ z (-^b)^ + (^b) 



1 



(C.13) 



where cu B = w^- — u Fl 



27r( ^ g ^ and we define JC ijU)B and X iA , B by 



and 



K 



l,UJ B 



det G^.det G LUFp+u)B 



1 d 

2J_ 



#■ tanh § <p,- 
— tanh ^(p.; 



, when ujb = 0, 
when c<jb 7^ 0. 



(C.14) 



£r tanh + ^- tanh |<^j , when cjb = 0, 
g|rtanh f ^ - (^) 2 K ijUlB , when ^ B ^ 0. 



(C.15) 



Gathering all terms of second order with respect to the fluctuations 5(p from ()4.43j) one 
gets (|4.52|) if the mean-field (p is of a Neel or ferromagnetic type (ip + = (p~ = and 



C.3 Derivation of the free energy with fluctuation 
contributions 



From relations (|4.52|) one can define 
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®= (ftanh'fen ) (C.16) 



4 V 2 



and 



^ /ltanh (§£>f)' 
y2 iff - iu B 

Using relations (|4.55J) and (|4.56J) with rhi = (— l) % fn + AmBi, ® and (2) read 



(C.17) 



® = (ftaah'gfc 
= f (1 - ) 



[la] + [16] (-1)*-* 

f (l - 4 (m + Am .B) 2 ) , when ^ = (-1)***-* 

f (1 - 4 (m 2 + (Am .5) 2 )) - 2(3Am x0 .Brh( y -lf^ , when 5, = B. 

(C.18) 



and 



= A tanh(fO 
\2 tp z a - iu B 

where 



r2l (m + Amo.B).(B + 2D| J|(m + Am .ff)) 

L JwB [(5 + 2D|J|(m + Am J)) 2 + w B 2 ] 1 J 

roftl = iu(m + Am .B) 

B [(B + 2D\J\{rh + Am .B)) 2 +u B 2 ] 1 ' 



if 5j = and 



Am x0 .B (B - 2\J\B.Am x0 - iuo B ) - 2\J\rh 2 



L luJB {B - 2\J\Am x0 .B) 2 - {2\J\mf V ; 

rh(B — 2\J\B.Am x0 - iu B ) - 2|J|m.Am x0 .fl 2 
[ Ls " (B - 2\J\Am x0 .B) 2 - {2\J\mf 1 ' ' 
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if B, L = B. After these developments in equations (|4.52|) . the use of a Fourier transfor- 



mation and integration over the Hubbard-Stratonovich auxiliary field S(p equations (J4.4 
and (|4.5()jl come out as 



Z 7 



n det 



keSBZ 



z +- = n n det 



w b keSBZ 



-1/2 

1 - [la] J{k)\ - [16] J(k + tt) 
- [16] J(jfe) (l - [la] J(k + tt) 

1 - [2al B J(k)) - [2bl B J(k + tt) 
- [2bl B J(k) (l - [2al B J(k + tt) 



-i -i 



(C.23) 



(C.24) 



Here J(k) = — Z\ J\ju and n is the Brioullin vector relative to the spin lattice as defined 
in section l4~2l Then free energy components Tmf-, Fzz and are given by 



,2 , jp 



J~mf = AfD\J\(m + Am .B) z -^-\ncosh[^-[B + 2D\J\(rh + Am .B)} 



(C.25) 



iV In 

2(3 ^ 



keSBZ 



I - I (3D ^ ) [l-A(m + Am .B) 2 ]- 



(C.26) 



r + - =| £ - 



sinh ([5 + 2D | J|(m + Am .B)] 2 - [2D\J\^(m + Arh .B)]' 
sinh (f [5 + 2D\J\(m + Am .B))) 

if the magnetic field Bi is equal to (— l) ri ' n B and 



1/2 



(C.27) 



T MF =HD\J\ (rh 2 - (Am x0 .B) 2 ) 



- ^ In | cosh 



- ((1 - 2D\ J\Am x0 ).B + 2D\J\m) 



. cosh 



^ ((1 - 2D\J\Am x0 ).B - 2D\J\m) 



(C.2* 
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T 

J 7. 



h E ln I 1 - (^Y^) 2 t 1 - 4 ^ - Am xo-B) 2 } [1 - 4(m + Am x0 .B) 2 ] 

(C.29) 



2(3 



keSBZ 



4 E >» 



sinh 2 (3D\J\fh< 1 -7 



1 - 



■ Am x0 .B x 2 



- sinh 2 f ((1 - 2D\ J\Am x0 ).B)- 



sinh 2 /3D\J\m- sinh 2 |((1 - 2D| J| Am x0 ).5) 



(C.30) 



if the magnetic field 5j is unform and equal to B. Here A/" is the number of lattice 
sites (sublattices A and B). The magnetization of the spin system is obtained by means 
of (jC.25|) . (jC.26|) and (jC.27|) . The magnetic spin susceptibility can be worked out with 

d2ni) and ra. 



C.4 The free energy of the XXZ-model with a stag- 
gered magnetic field 

The free energy of the XXZ-model in a staggered magnetic field Bi = (—l) l B is given by 
the three contributions Tmf-, T zz and where 



Tmf = ND\J(l + 5)\(m + Am xxz .B) 2 

- j lncosh (J [B + 2D\J{1 + 6)\(m+ Am xxz .B)]j 



(CM] 



zz = 2/3 2^ 1, 9 ) L 1 ~ 4 ( m + Am xxz.5) j 



keSBZ 



(C.32) 
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=| £ In 



sinh (§ ([5 + 2,01/(1 + 5)1(771 + Am^.B)] 2 - [2£>| J\ 1% {m + Am xxz .£?)] 2 ) 



1/2 



sinh (f [5 + 2£>| J(l + (J)|(m + 

(C.33) 



Appendix D 

Diagonalization of Mean-Field 
Hamiltonians 



D.l Bogoliubov transformation on the Neel mean- 
field Hamiltonian 

The Neel mean-field Hamiltonian (|5.11|) 



n 



MF 



id 

+ E(4t fl 



-/i + f 



2 



ft 
2 



fi,1 
fi,l 



(D.l) 



can be diagonalized by means of a Bogoliubov tranformation Introduce the linar 

combination of new creation and annihilation fermion operators /3// ± ^ and A,(±) 



where the unitary matrix U is define by 



At \- n ( &<+) j (D-2) 



17 



it, 



(D.3) 



Coefficients {-Uj} and {t>j} have to verify det U = 1 and {7 1 = W in such a way that 
in the transformation ()D.2|) the creation and annihilation operator {0l(±\} and {/3j,(±)} 
anticommute. To diagonalize the Hamiltonian (|5.1()j) one should find the coefficients {u^ 
and {v^ for which 
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U^HmfU = U ] H MF U 



h3 



t 

«,(-) 



By means of the transformation 



(D.4) 



Mi = e loli cos 0j 
Ui = sin 9i 

where and ^ are angles, WHmfU is diagonal for 



(D.5) 



tan a,- 



tan 29j 



-y 



\lm 2 + m 2 



(D.6) 
(D.7) 



And then 



U^HmfU 



uj 



(PFP) 
*,(+) 






(PFP) 
*,(-) 



(D.8) 



with a; 



(PFP) 



/i ± (PFP refers to the Popov and Fedotov procedure). 



D.2 Bogoliubov transformation on Diffuson mean- 
field ansatz 

The mean-field Hamiltonian ()5.39|) with the imaginary chemical potential fixing exactly 
the number of spin per lattice site given by 



n 



(PFP) 
MF 



*>m + EE (4/, 



keSBZ 



-fi + A cos \z^ kx ,k v -«Asin|^7 fea!i A 



-MA sin^7 fe 



-ft — A cos f 27^ 



f-. _ 



(D.9) 
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can be diagonalized by the transformation (jD.3 



J k+1T,tT 



ut v; 



The transformation matrix 



U 



is unitary and 



(D.10) 



(D.ll) 



KH 2 + KH 2 = i 



(D.12) 



so that /3* and (3 are fermion creation and annihilation operator. Following the section 
ID. II we define 



Ut 



e^i cos 9- k 



The mean-field Hamiltonian is diagonal if 

tan 2% 

where 7^ = | (cos k x + cos k y ) and reads 



7T 

2 

-fey 



(D.13) 
(D.14) 



(D.15) 
(D.16) 



zA 5 

T7T 



+ EE[ 



; (PFP) /? f 

(+),fe,<T M (+),fc,<T P ( + )^>' 



keSBZ 



with 



(D.17) 



; P f r P) = -fi + 2AJcos 2 k x + cos 2 k 

(PFP) 
(-),£,er 



-/i — 2 A a/ cos 2 /c z + cos 2 



(D.18) 
(D.19) 
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D.3 Bogoliubov transformation on Cooperon mean- 
field ansatz 



The mean-field Hamiltonian (|5.fj()j) describing the cooperon mean-field ansatz with exact 
occupation parameter /i is given by 



and rewritten as 



keBZ,a 



k,(T —k, 



n 



(PFP) 
MF 



AT 



\J\ 



keBZ 



where 



^ll,k,a ^12,k,a 
^2l,k,a ^22,k,a 



h 



k.—a 



;d.2o) 



(D.21) 





f 

2 


H-12,k,a = 


zV lk 
-a - 

2 


^ 21, k,tr = 


zV lk 
—a - 

2 


^22,k,a = 


A* 
2 



(D.22) 

(D.23) 

(D.24) 
(D.25) 



Following the same step as in section ID .21 one define the Bogoliubov transformation 



A;, cr 



/ 1 



k. — a 







k,—a 



;d.26) 



where [/ is the unitary matrix 



U 



u k v k 



The coefficients verify \uA 2 + \vA 2 = 1 and 



(D.27) 



e^cos% 
sin % 



(D.28) 
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In order to diagonalize the mean-field Hamiltonian one looks for angles 6^ and p^ such 
that 



u-^Pu = ( PL A 



Developing U ^mf^U in detail one gets the system of equations 



k.a 



10 



(PFP) 
(+),k,a 

u) 





(PFP) 
(-),k,a J 




(D.29) 



r® 



H n utuj; - H X2 u%v % - H 2l up)t + H 22 v* t v 



b k v k 



k v k 



CO 



(PFP) 



k k ~(+) I jt ) < T 
(PFP) 



H llU tvt + H 12 (utf - H 2l {v^ - H 22 vtut = ' 



Introducing (jD.28|) in (3) and (3) leads to 



(D.30) 



H n e~ ilp K2 cos % sin % + F 12 (e~ 2 ^ cos 2 % - sin 2 %) = (D.31) 
Hne^.2 cos % sin % + F 12 (e 2i ^ cos 2 % - sin 2 %) = (D.32) 



Separating real and imaginary part one gets 



Tie part Xm part 

H 12 cos 29^ cos <£>g = and Xm ^ ) sin 2#g — Hi 2 sin = 
i?i2 cos 26 j: cos y?£ = and Xm ^' sin 2#£ + i/12 sin ip^ = 



(D.33) 



Equation (3) and (§) cannot be verified simultaneously. As a consequence there is no 
unitary matrix U and hence no Bogoliubov transformation which can diagonalize the 
mean-field Hamiltonian (|5.fi()|) . The mean-field Hamiltonian of the cooperon does not 



conserve the number of particles and hence gets in conflict with the Popov and Fedotov 
procedure which fixes the number of particles strictly 
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Appendix E 

Derivation of the QED% action and 
the polarization function at finite 
temperature 



E.l Derivation of the Euclidean QED action in (2+1) 
dimensions 

At low energy near the two independent points k — (±f , f ) + k of the Spin Brillouin 
Zone (see figure 14. 2|) the Hamiltonian ()6.5|) can be rewritten in the form 



H 



y y(f^ / f - 

r 3 



keSBZ a 
-/ill + V2A 



kfr. 



r 3 



kyl 



+ V2A 



-k r I q _° r2 ) +/>',. /B 



/ fl,k,a \ 
f l,fc+7r,<r 

f"2,%,(T 

\ f2,k+n,a / 



(E.l) 



with 7T = (7r, 71") the Brillouin vector. r 1; r 2 and r 3 are Pauli matrices 



-1 :;).*-(! j «-» 



f ' - and f , r ( f - and f„ ? ) are fermion creation and annihilation operators near the 

J l,k,a J l,k,cr \->2,k,a J 2,k,cr/ r 

point (f,|) ((-f,f)). 
Rotating the operators 
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fk ~ y/2 ija$ + fb,k 
fk+ir ~ "73 (ja,k ~ fb, 



leads to 

"= EDI 

keSBZ a 

where k+ = k x + k y and k. 



fil + Ak^ 



Tl 



r 2 

k x — k y , A = 2A cos | and 



Ak_ 



r 2 

Ti 



( f la,kcr \ 
f lb,kcr 
f, 2aka 
\ $2bko I 



In the Euclidean metric the action reads 



S h 



+ iAk 



keSBZ a 
T 2 

+ 1 



Tl 



r 3 
r 3 

iAk_ 



Tl 

-r 2 



r 3 
r 3 



Through the unitary transformation 



ka 



1 

e^ T3 



1 

-Ti 



ka 



and writing k + = k 2 and A;_ = ki 



keSBZ a 

where = ip^j and the 7 matrices are defined as 



7 



r 3 



, 7 







, 7 



r 2 
-r 2 



(E3) 



(E.4) 



(E5) 



(E6) 



(E.7) 



(E8) 



(E9) 



-t 3 ; ' ' v -n 

Using the inverse Fourier transform ijj^ = J d 2 rip? cr e lk ' r the Euclidean action reads finally 

p/9 



5 B = / dr / rf 2 r^^ 7 («9 T - //) + A 7 fc <9 1 



(E.10) 



With a "light velocity" = (1, A, A) leading to a curved metric as explained in subsection 

EM 
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E.2 Derivation of the photon polarization function at 
finite temperature 

The Fourier transformation of the second term of the spinon action given by equation 
(j6.15J) reads 



sf[i>A 



EE 



(2tt) 2 /3 



(2vr) 2 

Hh-h)- 



d 2 k 2 T 



i>a (fcl 



(2vr) 2 
(2n) 2 f3) 2 



(E.ll) 



with k = (up 



2n , 



n + 1/4), A;). Integrating over the fermion field ip and keeping the 



second order in the gauge field leads to the effective gauge action 



Sf] f [o] = ~Tr [Gp.igYaf (E.12) 



2 



with Tr = I j0^2-^2 I 7^f^ r - The trace tr extends over the 7 matrix space, and 
G F 1 (ki — k 2 ) = ij2^^fi(ki — k 2 ). The pure gauge action comes as 



1 \r^\ ^ r d 2 ki 1^ 



(2«rp^j (2tt) 



^.7%(^-^").5.7^ 



jfe" 



With the change of variables k\ — k = q and k\ = k 



S 



(2) 
eff 



-y 



d 2 q 
(2k) 



2 «7M 



-g)n^(g)a,(g). 



where q = (lub = ^S-rn, q) and the polarization function is given by 



n-w - J EE / 



a top 



k 2 (k + q) 



(E13) 



(Ell) 



(E15) 



Then using the Feynmann identity i = dx ( ax+( ^_ x ^2 ^ v can be rewritten as 
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By means of a change of variables k — > k' — xq and using the identity tr [yy^-y 1 '] 
4. [5pn.8 vv - 5 pr) .5^ + 5 p „.8 m ] one obtains 
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where a = g 2 ^]l. Following Dorey and Mavromatos Lee |1H], Aitchison et al. [Sj 



(T=l 



and Gradshteyn [SO] we define 



Si 



So 



s* 



E 



1 



n =-oo W + z(l - 
f3 2 [ sinh(27rF) 



E 



cosh(27rE) - cos(2ttX) 
1 



(3 2 1 dS x 



+ x(l — x)q 2 



E 



a;, 



8tt 2 ' Y dY 
(3 dSx 



oo \ k' + x(l — x)q 2 



4vr<9X 



(E.18) 
(E.19) 

(E.20) 



with X = x.m + 1/4 and Y = -^y k' + x(l — x)q 2 . The polarization can be expressed 
in terms of these sums and reads 
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for the spatial components. 

— * 

Integrating over the fermion momentum k' one gets 



where 



n 00 = n 3 -4ni-n 2 (E.23) 
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n 3 = — / rfxln2 J D(X,r) (E.27) 



Jo 

and D(X, Y) = cosh ^/3g-\/:r(l — x)j + sin(27ra;m). 
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